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^ ! We consider stability of regimes of hydromagnetic thermal convection in a rotating 
horizontal layer with free electrically conducting boundaries, to perturbations involv- 

Q ' ing large spatial and temporal scales. Equations governing the evolution of weakly 

' nonlinear mean perturbations are derived under the assumption that the a— effect 

■ is insignificant in the leading order (e.g., due to a symmetry of the system). The 

^ ! mean-field equations generalise the standard equations of hydromagnetic convection: 
New terms emerge - a second-order linear operator representing the combined eddy 

^ \ diffusivity, and quadratic terms associated with the eddy advection. If the perturbed 

J> ' CHM regime is non-steady and insignificance of the a— effect in the system does not 

^ ■ rely on the presence of a spatial symmetry, the combined eddy diffusivity operator 

! also involves a non-local pseudodifferential operator. If the perturbed CHM state is 

^ \ almost symmetric, a— effect terms appear in the mean-field equations as well. Near 

'nI" I a point of a symmetry-breaking bifurcation, cubic nonlinearity emerges in the equa- 

^ ' tions. All the new terms are in general anisotropic. A method for evaluation of 

O ■ their coefficients is presented; it requires solution of a significantly smaller number of 

>■ . auxiliary problems than in a straightforward approach. 

^ ' 1. Introduction 

An attribute of convection in a melted planetary core is the presence of interacting 
structures involving different spatial scales. One example of such a structure is the 
Ekman layer in rotating convective flows, whose instability can result in magnetic field 
generation (see Ponty et al, 2001a,b, 2003; Rotvig and Jones, 2002). Core-mantle 
coupling, which is supposed to be responsible for the length of day variation at time 
scales of decades, is another one: e.g., topographic coupling is due to topographic 
structures at the core-mantle boundary, which do not exceed 5 km in size (Bowin, 
1986; Merrill et al., 1996); this is small compared to the radius of the liquid core (see 
also a discussion of the influence of irregularities of the boundary of the Earth's liquid 
core on the flow in it and on the Earth's magnetic field in Anufriev and Braginsky, 
1975, 1977a,b). While in these two examples the small-scale features are located at 
the boundaries of the liquid core, scale separation can occur in the entire volume, 
where convection takes place. For instance, in geostrophic flows in rapidly rotating 



1 



spherical or cylindrical shells fluid moves in the so-called Taylor columns, parallel to 

the axis of rotation; the columns have a smaller width than the size of the container 
of the fluid. Narrow (in the horizontal direction) cells emerge in thermal convection 
of fluids rapidly rotating about the vertical axis (Bassom and Zhang, 1994; Julien 
and Knobloch, 1997, 1998, 1999; Julien et al, 1998) and in magnetoconvection in 
the presence of strong magnetic fields (in the limit of large Chandrasekhar numbers; 
see Matthews, 1999; Juhen et al, 1999, 2000, 2003). A hierarchy of scales and 
related phenomena, such as intermittency and direct and inverse energy cascades, 
are characteristic features of turbulence (Frisch, 1995). 

This suggests that analytical asymptotic methods can be used in combination with 
numerical ones to study various convective and MHD regimes. This approach was 
followed, in particular, in the well-known geodynamo studies by Braginsky (1964a-d, 
1967, 1975) and Soward (1972, 1974), who constructed asymptotic expansions of 
solutions to dynamo problems to determine the values of the magnetic a-effect coef- 
ficients. 

In multiscale MHD stabihty problems it is typically assumed that the characteris- 
tic spatial and temporal scales of a perturbation are much larger than the character- 
istic scales of the perturbed regime. We assume that the perturbation depends on the 
so-called fast spatial x and temporal t variables, and on the slow ones^: X = s{xi, X2), 
T — sH. The small parameter e is the spatial scale ratio. In what follows, vector 
fields which depend only on small-scale (fast) variables will be called small-scale fields, 
and those depending as well on large-scale (slow) variables will be called large-scale 
fields. In this terminology, we are concerned with stability of small-scale regimes to 
large-scale perturbations. 

The perturbation is expanded in power series in the small parameter e. General 
methods of the theory of homogenisation for equations in partial derivatives for multi- 
scale systems (see, e.g., the monographs Bensoussan et a/., 1978; Olcinik et al, 1992; 
Jikov et al., 1994; Cioranescu and Donato, 1999) are then applied to derive rigorously 
a closed system of equations for large-scale structures in the perturbation, averaged 
over small scales. The influence of the small-scale dynamics is represented by new 
terms (often called eddy corrections) emerging in the equations. Computation of 
coefficients in these terms requires numerical solution of systems of linear differential 
equations in partial derivatives - the so-called auxiliary problems. The multiscale 
asymptotic techniques offer a tool to split analytically the large- and small-scale dy- 
namics in the complete problem, if the small-scale dynamics is in a certain sense 
homogeneous. Consequently, a full resolution is not required for the solution of the 
auxiliary problems, which is a significant advantage of the combined approach. The 
techniques were applied to specific two- (Sivashinsky and Yakhot, 1985; Sivashinsky 
and Frenkel, 1992; Novikov and Papanicolaou, 2001; Novikov, 2004) and general and 
specific three-dimensional equations of hydrodynamics (DubruUe and Frisch, 1991; 
Wirth, 1994, Wirth et al, 1995), passive scalar transport (Biferale et al, 1995; Ver- 
gassola and Avellaneda, 1997; Majda and Kramer, 1999) and kinematic magnetic 

^The width L of the layer is fixed, and hence no slow spatial variable in the vertical direction is 
introduced. 
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dynamos (Lanotte et al, 1999; Zheligovsky et al, 2001; Zheligovsky and Podvigina, 
2003; Zheligovsky, 2005). The effect of eddy viscosity was observed in direct numer- 
ical investigation (Murakami et at, 1995) of stability of planar flows to large-scale 
three-dimensional perturbations (without a prior recourse to asymptotic analysis). 

Nepomnyashchy (1976) considered eddy viscosity in the Kolmogorov flow and 
found that near the bifurcation where eddy viscosity passes through zero, the regime 
satisfies an equation with a cubic nonlinearity of the Cahn-Hilliard type, studied by 
She (1987). E and Shu (1993) studied numerically inverse cascade in Kolmogorov 
flow, using homogenised nonlinear "effective" equations. Mean-field equations for 
weakly nonlinear perturbations of two-dimensional hydrodynamic steady flows were 
considered in the terms of the stream function by Gama et al. (1994) and Prisch et 
al. (1996) (in the latter paper an additional term describing the /9-effect was intro- 
duced into the Navier-Stokes equation to take into account rotation of fluid, and the 
case of marginally negative eddy viscosity was assumed). Newell (1983) and Cross 
and Newell (1984) investigated equations for perturbations, analogous by the method 
of derivation to the mean-fleld equations under consideration; however, their deriva- 
tion was carried out for model equations describing approximately convective flows in 
a layer in the form of deformed rolls. Newell et al. (1990a,b, 1993, 1996) and Ponty et 
al. (1997) explored weakly nonlinear dynamics of convective patterns and the devel- 
oping defects in a system of rolls. They considered a complete system of equations of 
the Boussinesq convection in a layer with rigid boundaries, employing the variables 
the amplitude and the phase; the mean-field equation for the slow phase was derived. 
(No magnetic field was considered in the papers referred to in this paragraph.) 

Using the multiscale asymptotic techniques, Zheligovsky studied linear (Zheligov- 
sky, 2003) and weakly nonlinear (Zheligovsky, 2006a) stability of parity-invariant 
three-dimensional space-periodic MHD regimes, and Baptista et al. (2007) consid- 
ered linear stability of convective hydromagnetic (CHM) steady states in a layer, 
symmetric about a vertical axis. As originally found by Dubrulle and Prisch (1991), 
in this class of multiscale problems a term manifesting the a-effect appears in the 
mean-field equations in the leading order, s; no other terms, linear or nonlinear, sur- 
vive averaging at this stage (whichever kind of stability - linear or weakly nonlinear - 
is inspected). Thus, one obtains mean-field equations involving other "eddy" effects, 
which can give rise to non-trivial dynamics of mean fields, only if in the leading order 
s the a-effect is insignificant (e.g., if appropriate components of the a-tensor vanish). 
The symmetries of the perturbed regimes assumed in the papers mentioned above in 
this paragraph guarantee that the a-effect vanishes, or it is insignificant. However, 
the symmetries are not necessary for this: e.g., the AKA (kinematic a) effect does 
not emerge in ABC flows (Wirth et al, 1995). Zheligovsky (2006b) considered weakly 
nonlinear stability of three-dimensional CHM regimes in a layer just assuming that 
the a-effect is insigniflcant. 

In the linear stability problems mentioned above, the main terms of the series 
expansion of the perturbation modes of steady states and their growth rates are 
eigenvectors and eigenvalues, respectively, of the operator of combined eddy diffusiv- 
ity. This is a linear partial differential operator of the second order, which is not, 
in general, isotropic or negatively defined. If it has eigenvalues with a positive real 
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part, one says that there is negative eddy diffusivity (Starr, 1968). In the weakly 
nonUnear stabihty problem it is assumed that the amplitude of the perturbation is 
of the same order as the scale ratio, e, and the evolution of such perturbation due to 
unabridged nonlinear system of equations is considered. Then additional quadratic 
terms emerge in the mean-field equations for the perturbation, which are analogous 
to the advection terms of the original equations. 

Lanotte et al. (1999), Zheligovsky et al. (2001), Zheligovsky and Podvigina (2003) 
and Zheligovsky (2003, 2005) have shown that the presence of small scales is beneficial 
for the action of kinematic magnetic dynamos and for the growth of magnetohydro dy- 
namic instabilities, and the effect of eddy viscosity can provide a working mechanism 
for development of these phenomena. A question immediately arises: What happens 
to this mechanism, when the initial stage of development of an instability is over and 
nonlinear effects switch on? It is natural to search for preliminary answers to this 
question in the context of weakly nonlinear stability theory. 

In the present paper we are concerned with weakly nonlinear stability of regimes in 
Rayleigh-Benard convection in a layer in the presence of magnetic field. We consider 
the so called "forced" convection, where external forces (in addition to the buoyancy 
and Coriolis forces) are supposed to act, and/or sources of heat and magnetic field 
are present inside the layer (section 2). In their absence the system is translation 
invariant in horizontal directions, and as a result the problem becomes more complex; 
this CHM system will be considered in a sequel to the present paper. 

The leading term in the expansion of perturbations involving large scales is com- 
prised of amplitude-modulated neutral modes of the operator of linearisation of the 
equations, governing the behaviour of the CHM system: 

(vo,ho,^o) = EC7'=s„ 

k=l 

where s^(x, t) are the small-scale neutral modes, and C'^(X, T) are amplitudes, slowly 
varying in space and time. Our goal is to derive the equations, describing the evolu- 
tion of the amplitudes C^. The derivation rehes on existence of constant vector fields 
in the kernel of the operator, adjoint to the operator of linearisation near the per- 
turbed regime. The homogenised equations are solvability conditions for equations in 
fast variables, which by the Fredholm alternative theorem (see section 4) amounts to 
the orthogonahty of the right-hand sides of the equations to the kernel of the adjoint 
operator, and in our case is equivalent to averaging of the respective components of 
the equations. Direct substitution shows that any combination of a constant fiow 
and magnetic field with a zero temperature component belongs to the kernel, if they 
satisfy boundary conditions. Thus the boundary conditions determine the number of 
the homogenised equations, and the auxiliary problems to be solved in order to eval- 
uate the coefficients of the eddy terms in these equations. For electrically conducting 
free boundaries held at fixed temperatures, considered here, the dimension K of the 
kernels of the linearisation and of the adjoint operator is at least 4. 

In the problem under consideration, generically K = 4, and mean horizontal com- 
ponents of the flow velocity and magnetic fleld of the neutral small-scale modes 
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are non-zero. Hence we consider the neutral modes = (S^'" + e^, S^' , S^' ) and 

— i^k'^j^k''^ + ^fcjSfe'^), k = 1,2, where horizontal flow velocity and magnetic 
fleld components of all vector flelds S^' vanish. Let ((f)) denote the average over fast 
variables of the horizontal component of a vector field f . Clearly, 

(( E C% }^ = (((vo))., 0, 0), ( E C% }^ = (0, ((ho))., 0), 

fe=l fe=3 

i.e., are the mean horizontal components of perturbations of the fiow, ((vo))/i, and 
magnetic field, ((hg)).. Consequently, we call "mean- field" the homogenised equations 
for C*^, that we construct. (Mean vertical flow and magnetic fleld components of 
perturbation vanish due to solenoidality and the boundary conditions.) 

In general, the a-effect emerges in the system (section 6). Then the mean-fleld 
equations turn out to be linear, resulting in instability of the CHM regime to large- 
scale perturbations. It can be, however, insignificant in the leading order, for instance, 
due to the presence of certain spatial or spatio-temporal symmetries in the system 
(section 7). Assuming it to be insignificant, as in Zhcligovsky (2006b), we derive 
(sections 8 and 9) the mean-field equations (59) and (62) for the evolution of the 
mean part of a perturbation, which generalise the standard equations of hydromag- 
netic convection. The linear operator of combined eddy diffusivity and quadratic 
terms associated with eddy advection emerge in the mean-field equations, like in the 
studies cited above. If the perturbed CHM regime is non-steady and the symmetry 
responsible for the a-effect insignificance in the leading order is spatio-temporal, non- 
local operators representing new physical effects emerge in the mean-field equations. 
We show (section 10) that if the perturbed CHM fiow is asymptotically close to a 
symmetric one, and its antisymmetric part is of the order of the scale ratio e (in par- 
ticular, this happens when a branch of CHM regimes emerges in a symmetry breaking 
Hopf bifurcation), then a-effect terms also appear in the mean-field equations, which 
now take the form (71) and (73). 

Finally, in section 11 we consider a symmetric CHM system close to a point of 
a symmetry- breaking pitchfork bifurcation for overcriticality of the order of e^. In 
this case the kernel of the operator of linearisation is five-dimensional {K = 5); in 
addition to neutral small-scale modes with non-zero horizontal components of the 
fiow velocity and magnetic field, there exist a mean-free neutral small-scale mode, 
whose amplitude (denoted cq in section 11) is not associated with any mean field. 
The mean-field equations (86) and (88) now constitute a closed system together with 
the equation for the evolution of this amplitude (since this equation is bulky, it is 
moved to appendix F). Eddy diffusivity, eddy advection and a-effect terms and non- 
local operators can still be found in the mean-field equations. Cubic nonlinearity 
emerges in the equation for the amplitude of the neutral mean-free mode, and if the 
perturbed CHM regime is non-steady and the symmetry responsible for the a-effect 
insignificance in the leading order is spatio-temporal, the equation involves a variety 
of linear and nonlinear non-local operators. 

Pairs of the mean-field equations (59), (62); (71), (73); and (86), (88) combined 
with the amplitude equation presented in appendix F, arising in the analysis of weakly 
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nonlinear stability of a single CHM regime and branches of regimes emerging in 
the Hopf and pitchfork bifurcations, respectively, constitute the main results of the 
present paper. Following Zheligovsky (2006b), we do not assume that the perturbed 
CHM regime is periodic in horizontal directions (although this geometry is convenient 
for numerical solution of the auxiliary problems) . This condition must be relaxed in 
order to consider stabihty of semi-regular structures, such as spiral-defect chaos in 
thermal convection (see Figs. 5b,d,f in Bodenschatz et al, 2000). Also, like Zhe- 
ligovsky (2006a,b), we do not assume that the perturbed CHM regime is steady. 
Consequently, in derivation of the mean-field equations averaging must be performed 
over the entire domain of the fast (small-scale) spatial and temporal variables. The 
coefficients in the newly emerging terms are constants, like this is the case if a steady 
CHM state is perturbed. (In particular, temporal dependence of the eddy diffusiv- 
ity tensor, studied by Gama and Chaves, 2000, is consequential only on time scales, 
which are not involved in the derivation of the mean- field equations.) 

In appendix G a method for evaluation of the coefficients is presented, which 
requires to solve a significantly smaller number of auxiliary problems than in the 
straightforward approach. This method relies on solution of auxiliary problems for 
the adjoint operator, like in Zheligovsky (2005, 2006a,b). 

Like Zheligovsky (2006b), we consider convection of fiuid, rotating about a vertical 
axis (which is more astrophysically sound than the absence of rotation). This gives 
rise to the following algebraic difficulty. In the presence of the Goriolis force, the 
kernel of the adjoint operator includes constant vector fields with a non-zero constant 
flow only, if a linear growth in horizontal directions of the potential of the subtracted 
gradient is allowed. However, in this case surface integrals appear from the pressure 
gradient when averaging of the equations over fast variables is performed, which is 
technically inconvenient. To overcome this difficulty we consider (like Zheligovsky, 
2006b) the Navier-Stokes equation for vorticity. As a result, the equation for the 
perturbation of the flow emerges as the solvability condition for equations in fast 
variables at the order of s^, and not e"^, as usual. Moreover, solution of equations at 
the order of £" is preceded by derivation of the spatial mean {\'n)h from the spatial 
mean of the vorticity equation at the order of 

2. Equations of hydromagnetic thermal convection, boundeiry conditions 
and the Hneeirisation operator 

A CHM regime V, H, T, the weakly nonhnear stabihty of which we are considering, 
satisfies in the Boussinesq approximation the equations 



— = uV^n + V X (V X - H X (V X H)) + V X (V X re3 + pTe^ + F), (la) 



an 

'dt 



TyV'H + V X (V X H) + J, 



(lb) 



'dt 



kv'^t - (V • v)r + s, 



V-V = V-H = 0, 



V X V = 17. 



(Ic) 

(Id) 
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Here Y{x,t) denotes the velocity and f2(x, the vorticity of an electrically con- 
ducting fluid flow, H(x, t) magnetic field, T(x, t) temperature, t time, u, rj and k 
kinematic, magnetic and thermal molecular difi^usivities, respectively, r/2 angular 
velocity of the rotation, PTe^ the Archimedes (buoyancy) force, is the unit vector 
along the coordinate axis Xk, F(x, t) a body force, J(x, t) reflects the presence of ex- 
ternally induced currents in the layer and 5'(x, t) of heat sources. To simplify the no- 
tation wc will occasionally use 10-dimensional vectors (w, v, h, 9) and 7-dimensional 
vectors (a;,h, ^^) or (v,h, ^). 

The following conditions are assumed on the horizontal boundaries of the layer: 
• free boundaries: 



dVi 



X3=±L/2 



dx:\ 



= 0, 



xz=±L/2 



X3=±L/2 



= 



(2a) 



X3=±L/2 

electrically conducting boundaries: 



xz=±L/2 



dHi 



dx?, 



X3=±L/2 



dx?, 



X3=±L/2 



0, 



0, 



on. 



dxa 



0; 



X3=±L/2 



X3=±L/2 



0; 



(2b) 



(2c) 



boundaries at flxed temperatures: 



r 



X3=-L/2 



X3=L/2 



(Vector components are enumerated by the subscript.) It is convenient to introduce 
a new variable © = T — 7^ -|- b{x-i -\- L/2), where 5 — (7^ — '7^)/L (convection is only 
possible for 5 > 0). It satisfles the equation 



90 
'dt 



kv^q - (V • v)e + SV3 + S 



and homogeneous boundary conditions 

e 



X3=±L/2 



= 0. 



(3) 



(4) 



Linearisation of (1) in the vicinity of the CHM regime V, H, is the operator 
£ = {D'.C'.C'^), where 



C^{uj, V, h, e) = + i/V^w + Vx(Vxu; + vxl2 



dt 



-H X (V X h) - h X (V X H)) + r— + PV0 x eg, 

OXs 



dh 



jC^{v, h) = -— + rjV^h + Vx(vxH + Vxh), 

C/ C 



(5a) 
(5b) 
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BO 

C%v, 9) = -— + kVH - (V • V)^ - (v ■ V)0 + ^^^3, (5c) 

and a; = V X V. Evidently, = (Vx)£", where 

h, e,p) = -— + vV\ + Vx(Vxv) + vxJ2 
at 

-H X (V X h) - h X (V X H) + TV X 63 + pOez - Vp 
is the hnearisation of the Navier-Stokes equation for the flow velocity: 

= i/V^V + VxO-Hx(VxH)+rVxe3 + (5Qe^ - VP. 

at 

Here P = P' / p + |Vp/2 + T^(xf + x|)/8 is a modified pressure, P' pressure and p 
fluid density. 

We assume that in a weakly nonlinear regime the amplitude of a perturbation is of 
the order of e. The perturbed state fi + eo;, V+£v, H + eh, Q + ed satisfies equations 
(l)-(4), and hence the profiles of perturbations a;,v,h, ^ (which we henceforth call 
perturbations) satisfy the equations 

C^iuj, v, h, ^) + £V X (v X w - h X (V X h)) = 0, (6a) 

d'{w, h) + ffV X (v X h) = 0, (6b) 

£^(v,^)-£(v V)^ = 0, (6c) 

V ■ h = 0, (6d) 

V ■ a; = V ■ V = 0, (6e) 

V X V = a;. (6f) 

Define the spatial mean over fast variables and fluctuating parts of a scalar or 
vector field / as 

1 fL/2 ,112 re/2 

if) = lim — / / / /(x) dx, dx, dxs, {/} = /- (/), 
£-^■00 J-L/2 J-e/2 J-e/2 

and the spatio-temporal mean and the fluctuating part of / as 

((/)) ^ hm - r(/(x, t)) dt, ifi^f- ((/)), 

t^oo t Jo 

respectively. Denote by {f)k and {f}k the A;— th components of the means (f) and 
((f)), respectively, 

(f), = (7)363, {f}, = f-(f),; 

((f)).^((/)>3e3, OT.^f-((f)).; 

(f),= (/)l6i + (/)2 62, {f}, = f-(f),; 

((f)), ^ ifhe, + ifhe,, mn ^ f - ((f)).- 

(Thus the subscripts v and h are used to denote vertical and horizontal components 
of three-dimensional vector fields; we will also use the superscripts v and h to denote 
the flow velocity and magnetic field components of 7- or 10-dimensional vector fields.) 
Averaging the horizontal component of the equation for the fiow perturbation. 
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h, e,p) + £(v X (V X v) - h X (V X h)) = 0, 



over the layer of the fluid, find 

^ = (v), X res - (Vp).. 

Thus, the mean horizontal component of the perturbation of the flow is controlled by 
the mean of the gradient of the pressure perturbation p, which must be specified. We 
assume that there is no pumping of fluid through the layer due to pressure variation 
at the infinity, i.e. the growth of pressure at the infinity is bounded so that {'Vp)h = 0. 
Under this condition, if initially the horizontal component of the perturbation of the 
velocity vanishes, then at any time 

(v). = 0. (6g) 



3. Asymptotic expansions of large-scale weakly nonlinear perturbations 
of CHM regimes 

We introduce the slow spatial, X = e{xi,X2), and temporal, T = e'^t, variables. The 
exponent in the temporal scale ratio is tailored for CHM regimes, where the a— effect 
is insignificant in the leading order. A solution to the problem (6) is sought in the 
form of power series: 

oo 

u;=5:u;„(x,i,X,T)£", (7a) 

n=0 
oo 

v = ^v„(x,i,X,T)£", (7b) 

n=0 
oo 

h=^h„(x,i,X,r)£", (7c) 

n=0 
oo 

e = ^^„(x,i,x,r)£". (7d) 

n=0 

Consider the series in e resulting from substitution of (7b) and (7c) into (6d) and 
(6e). The mean and fluctuating parts of equations at order n > are 

Vx • (v„)ft = Vx • {K)h = 0, (8a) 
Vx-K},, + Vx-{v„-i},. = 0, (8b) 
Vx ■ {K}h + Vx ■ {K-i}h = 0. (8c) 

(In differential operators with the indices x and X differentiation in the respective 
spatial fast and slow variables is performed; Vx = {d/dXi,d/ 8X2,0). Henceforth 
differentiation in fast variables is assumed in the definition (5) of the operator C. All 
quantities for n < are zero by definition.) Clearly, 

Vx • {uJn)v = 0, (8d) 
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and V ■ u) — implies 

Vx • Mv + Vx • {oJn-i}v = (8e) 

for all n > 0. 

Evidently, {uJn)v, {'^n)h and {hn)h satisfy the boundary conditions (2). Due to 
(8a), we can introduce a stream function X2, t, T) for {vn)h- 

Substitution of the series (7a) and (7b) into (6f) yields 

Vx X {Vn}h = - Vx X v„_i. (9) 

Consider this as an equation in {Vn}h- It can be shown using (9) after the change 
of index n n — 1 and (Be), that the right-hand side of (9) is solenoidal in fast 
variables. Consider the operator curl which acts from the space of solenoidal vector 
fields satisfying (2a) into the space of solenoidal vector fields. The domain of the 
adjoint operator, also a curl, is the space of solenoidal vector fields satisfying (2b); 
its kernel consists of constant vector fields (0, 0, C). Thus (9) has a solution as long 
as the average of the vertical component of the right-hand side vanishes^: 

{u;n)v = Vx X (V„_i)ft. (10) 

(For n = this reduces to {uJo)v = 0.) In the terms of the stream function this 
equation is V^ipn-i = {<^n)3- One can uniquely determine {^n-i)h from this equation 
under the condition that ■0n-i is globally bounded (so that (6g) is satisfied), if the 
mean of (a;„)^ over the plane of slow variables vanishes. 
In view of (10), (9) reduces to 

Vx X {v„}fe = {a;„}^ - Vx X {v„_i}fe. (11) 

After a substitution {^n}h = v -|- Vx-S, where S is a globally bounded solution to 
the Neumann problem 



dB 

VlB = -Vx • {Vn-l}h, 



= 0, 

X3=±L/2 



the system of equations (11) and (8b) takes the form of (6e)-(6g). 

A solution to the system (6e)-(6g) is v = V{A.}h, where A solves the Poisson's 
equation 



V^A = -V X u;, ^ 

ox?. 



X3=±L/2 



= 0, As\.s=±L/2 = 0, 

X3=±L/2 



^This condition is sufficient for a space-periodic in horizontal directions CHM system. In general 
a sufficiently fast decay near zero of the spectrum of the field, to which the inverse curl is applied. 
We assume that the examined CHM state is such that the inverse curl can be applied to any field, 
whenever the space average of the vertical component of the field vanishes. 
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and V is the projection of a three-dimensional vector field into the subspace of 
solenoidal fields: VA. = A — Va, a being a solution to the Neumann problem 

da 



= V ■ A, 



dx?. 



= 0. 

X3=±L/2 



This defines the operator 7?., inverse to the curl, which acts from the space of solenoidal 
fields satisfying {uj)^ = and boundary conditions for vorticity, into the space of 
solenoidal fields satisfying (6g) and boundary conditions for flows. 

We substitute (7) into (6a)-(6c) and transform them into power series in e. The 
equations emerging at different orders constitute a hierarchy (see (A1)-(A3) in 
appendix A); they can be solved together with (8) and (9) successively, by considering 
separately their mean and fluctuating parts. When considering the equations at order 
n, the mean vorticity equation at order n + 1 is also used to flnd {(jJn+i)v and (v„)/j. A 
closed system of nonlinear equations for the averaged leading terms of (7) (which we 
call the mean-field equations) will be derived as solvability conditions for equations 
in fast variables at the orders n — 2 and 3. 

4. Solvability of auxiliciry problems 

Vector fields Jl, V,H, O constituting the CHM regime, whose stability is examined, 
depend only on fast variables. We assume that they are smooth and globally bounded 
together with their derivatives^ and all average quantities in fast variables, considered 
below, are correctly defined. This condition is satisfied, for instance, if the perturbed 
CHM regime is periodic in horizontal directions and time (since then the domain, in 
which all the fields are defined, is compact). 

Auxiliary problems of the following structure will be considered: 

£-(a;,v,h,^^) = f-, £^(v,h) = f^ C%^,,e) = f, (12a) 

Vx-w = c^^, Vx-v = rf^ V^-h = d^, (12b) 

oji — Vx X V = a, (12c) 

and the means i^>j}v, ((h)) /i and {v)h are specified. Here the right-hand sides are known 
and the compatibility relations, obtained by combining (12b) with the divergences of 
(12c) and the first two equations in (12a), must be satisfied: 

Vx • a - d'^ = 0, (12d) 
+ yVlr + r^^V^-i^, (12e) 



+ 7?V2 d'* = Vx • f ^ (12f) 



dt 

a satisfies the boundary conditions for vorticity, 

(d-) = {d^) = 0. (12g) 



^This is understood in the following sense: a field f(x, X, T) is globally bounded together 
with its derivatives, if the field and all required partial derivatives are bounded, with the bounds 
depending only on the order of the derivative, slow variables and fast time. 
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Due to (12e) and (12f), it suffices that vorticity and magnetic field satisfy (12b) 
initially. 

The fields a;,v,h and 9 must be globally bounded with their derivatives and 
satisfy the same boundary conditions as r2,V,H and 0, respectively For vector 
fields from this class, the identities 

(£'^(a;,v,h,^)), = -^, (13a) 

(£V,h)). = -^, (13b) 

hold true, and hence 

((£-(a;,v,h,^))), = 0, {C^ iy ,h))) ^ = 0. (14) 

(14) implies that the conditions 

((P(x,i))), = 0, (15a) 

((f''(x,i))), = (15b) 

are necessary for existence of a solution to equations (12). Integrating (13) in fast 
time, find 

(a;),[^^ - = J\nvdt, (h),[^^ - (h), = J\f')hdt, (16) 

implying 

(a,), = (M), - {{ J^inv dt }} , (h), = ((hK - {{ J^{t')H dt } (17) 

(although the integrals in the right-hand sides fiuctuate only in time, the notation 
is still applicable). Thus, (a;)^ and (h)^ are well-defined, if and only if the 
means ((/q (f'^)t,(it)) and ((/q (f^)/i(it)) are. The latter conditions are stronger than 

(15) , since if ((f^))^ and ((f^))/i exist, (15) follows from existence of ((/q and 
iIo{^'^)hdt}. Averaging the vertical component of (12c) over fast spatial variables, 
find its solvability condition (see section 3): 

{^)^ = iu;h-{l\nvdt}. (18) 

In what follows we assume that for arbitrary globally bounded together with 
their derivatives smooth solenoidal zero-mean (15) fields f'^(x, t), f^(x, t) and /^(x, t) 
the system (12) has a solution u:,\,h.,6 in the considered class at least for some 
smooth initial conditions, satisfying (17) for t = 0. This is necessary to ensure that 
the auxiliary problems, stated in sections 5 and 8, have solutions with the required 
properties. A smooth solution to (12) can be constructed as a solution to a parabolic 
equation for smooth initial conditions and right-hand sides in (13a). However, it is 
not guaranteed that it is globally bounded, since the region occupied by the fiuid is 
not compact. 
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By virtue of (13), the spatial means (a;)„ and {h)h of solutions from this class to 
the system 

£(a;,v,h,^) = 0, (19) 

complemented by (6d)-(6g), are time- independent. A CHM regime is linearly stable 
to small-scale perturbations, if a solution to (19), (6d)-(6g) exponentially decays in 
time for any smooth globally bounded initial conditions satisfying (w)^ = (h)/j = 0. 
It is shown in appendix B, that if a CHM state is hnearly stable to small-scale 
perturbations, then for any smooth globally bounded right-hand sides in (13a) and 
initial conditions satisfying (12b), (12c) and (17), the solution to (12) is globally 
bounded. However, we do not demand that the perturbed CHM regime V, H, is 
linearly stable to small-scale perturbations (despite any small-scale instability has a 
larger, order zero, growth rate, than a large-scale instability), so that the formalism 
which is developed here were applicable to analyse weakly nonlinear instability of 
chaotic CHM attractors. 

In the remaining part of this section we show that, for a generic space-periodic 
steady or time-periodic CHM regimes, (15) are sufficient conditions for existence of 
a solution to (12), steady and/or having the same periods as the perturbed CHM 
regime. Define the operator of linearisation not involving the fiow velocity explicitly: 

M"^{uj', h', 6) = C{uj', 7^a;', h', 6), (20a) 

M"'{uj\ h!) = d'{nuj\ h'), (20b) 

M'\u:\e) = C\Tiw\e). (20c) 

The operator M' — (A4"^, M'^, A4'^) acts in the space of 7-dimensional vector fields 
{u}',h.',9), where vorticity and magnetic field are solenoidal, {(jj')v = 0, and the 
boundary conditions of the kind of (2b), (2c) and (4) are satisfied. Let Ai denote 
a restriction of M.' on the subspace in the domain of M. defined by the condition 
(h')/i = 0. Substituting 

w = + a, V = v' + V,>1^ + {Y)h, 

h = h' + V^A'^ + ((hK - { J\i'')h dt }} , 
where A" and A'^ are solutions to the Neumann problems 

dA'' 



dx-< 



= 0; ViA^ = d^, 



OA'' 



X3=±L/2 



0, 



X3=±L/2 



transform the problem (12) into an equivalent one: 

M^iu}', h', e) = f"*, 

MV,h',^) = r, 
Vx • w' = Vx ■ h' = 0, 



(21a) 
(21b) 

(21c) 
(21d) 
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(w'K^O, (v')/. = 0, (h')fe = 0, (21e) 
Vx • f"^ = Vx • f"" = 0, (21f) 
(f-), = (f"'), = 0. (21g) 

The operator £* = (C*)^, {C*Y) adjoint to £ = C", can be derived 

as usual by integration by parts in the defining identity 

(£*(v, h, 9) ■ (v', h', ^0) = ((v, h, 9) ■ C{V, h', 9')) : 

{c*y{w, h, ^) = — + i/vV - V X (V X v) 

+ P{H X (V X h) - V X ri - TV X 63 + 5^63 - 9VQ}h, 

dh. 

{C*f{w, h) = — + T^V^h + V X (H X v) + P(v X (V X H) - V X (V X h)), 

~ f)9 

(v, ^) = ^ + /^V^^ + (V • V)^ + Pv,. 

Conditions on horizontal boundaries for vector fields in the domain^ of C* can be 
obtained demanding as usual that the boundary surface integrals arising in this inte- 
gration are zero. For (2) and (4) defining the domains of C and C (together with the 
solenoidality condition for the flow, vorticity and magnetic components) , the domain 
of C* and C is the same. 

The operator M'* = {{M'*)^, {M'*)'\ (M'*)^), adjoint to M', can be found as 
an adjoint to a composition of C with the curl: 

(M'*r(oj, h,9)^^ + uV^uj + 7^'{-V X (V X (V X w)) 
ot 

+H X (V X h) + O X (V X uj) - r— + ^^63 - ^VG}/,, 

8X2. 

(A4'*)^(w,h) = — +r/V^h+Vx(Hx(Vxa;))+P((Vxw)x(VxH)-Vx(Vxh)), 

39 

{M'*)\u:, h,9)^ — + kVH + (V • V)^ + /3e3 • (V x w). 

Here the operator TZ' is also "an inverse curl" , like TZ, but it is defined for different 
boundary conditions: If w is a solution to the problem 

VxuJ^v-Vp ^ V^uj = -V X V, 

V • u; = 0, (uj)^ = 0, 

satisfying (2b), then TZ'v — uj. Since Al is a composition of M' with the projec- 
tion onto the subspace of fields with a zero spatial mean of the horizontal magnetic 
component, 

M*{uj, h, 9) = {{M'*r{'^. h, 9), {{M*)\<^, h, 9)}u, (M'*)^u;, h, 9)). (22) 
{M'* and M* have the same domains, as Ai' and M., respectively.) 

^Vector fields in the domains of C and C are assumed to be globally bounded with their deriva- 
tives; in this class the bilinear form (a ■ b) is not a genuine scalar product, since, for instance, 
(|ap) = for any smooth field a with a compact support. 
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Clearly, (0, (Ci, C2, 0), 0) e ker Al'* for any constant Ci and C2. We consider the 
generic case, where the kernel of A4'* consists of such constant vectors. (The case, 
where no external source terms F, J and S are present in (1), is not generic.) Then 
kevM* is trivial: for (u;,h, 6') G kerA^* 

M*(a;,h,^) = (0,(ci(t),C2(t),0),0) ^ M* (u:,h - J^J^ Ck{t') dt' e,,, 6^ = 0, 

by virtue of (22). Hence the property of ker AC* implies 

2 .t 

a; = 0, h-J2 Ck{i) diey, = (Ci, C2, 0), ^ = 0, 

and thus h = 0, since vector fields from the domain of Ai and M* satisfy (h)/j = 0. 

When the kernel of M.'* is two-dimensional, the conditions (15) are sufficient 
for existence of a steady or time-periodic solution to (21), which can be shown as in 
Zheligovsky (2003). This follows from the theorem on Fredholm alternative (see, e.g., 
Liusternik and Sobolev, 1961), stating that a solution to a problem {I+IC)X = ^ in a 
Hilbert space exists if and only if is orthogonal to ker(2r-|-/C*) . Here I is the identity 
operator, /C is compact and /C* is its adjoint. Apply the operators (—d/dt + z/V^)~^, 
{-d/dt + r/V^)"^ and {-d/dt + kV^)'^ to (21a)-(21c), respectively (By virtue of 
(21g) and the assumed boundary conditions they can be applied to the right-hand 
sides of (21a) and (21b); d/dt = in a problem with the steady data.) The system 
(21) is reduced thereby to an equivalent problem of the form M°{u}' , h', 9) — f", where 
M° acts in the domain of M. and has a trivial kernel, as long as ker A( is trivial. 
This problem is of the type considered in the Fredholm theorem. Thus application 
of the theorem guarantees existence of a unique solution with the required steadiness 
or time-periodicity. 

In the absence of space periodicity the condition (15) may be insufficient for 
existence of a globally bounded solution to (12), equivalent to (21). For instance, 
this can occur for CHM regimes, quasi-periodic in a horizontal direction: In this 
case, the Laplacian V^, acting in the space of vector fields the appropriate averages 
of which vanish, does not have a bounded inverse, and it is impossible to make a 
reduction of the problem, for which the Fredholm theorem can be readily applied. 

In sections 5-10 we assume that the problem (21) has a bounded solution for 
any right-hand sides f'^, f ^, f'^ (not imposing any periodicity conditions). Then the 
system of mean-field equations, that we derive, is comprised of equations for mean 
horizontal components of perturbations of the flow and magnetic fleld. In section 
11 we consider the case, where M* has an one-dimensional kernel; then the system 
of mean-fleld equations is expanded by a scalar equation for the amplitude of the 
respective mean-free neutral mode. 

5. Order equations 

The leading terms of the series (6a) -(6c) yield the equation 

/:(a;o,vo,ho,^o) = 0. (23) 
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Averaging over fast spatial variables the vertical vorticity component of (23), the 
horizontal magnetic component of (23) and the vertical component of (Al) for n — 1 
obtain, respectively, 

-^Qf- = (Ujo)^ = 

(provided {(jJo)v\t=o = to agree with (10) for n — 0), 

^^ = ^ (ho). = ((ho))., (24a) 

^^ = (c^i), = ((u;i))„. (24b) 

The last equation together with (8a) for n = and (10) for n = 1 implies 

(vo). = ((vo)). + v;(t,T). 

It is shown in appendix C that ( {vo}.|x=£(a;i,i;2)) is asymptotically smaller than any 
power of e. Thus (6g) asymptotically holds true to any order of e'" at any time, if 
v'q = 0, i.e. 

(vq), = ((voK, (24c) 

and the average of ((vq)). over the plane of slow spatial variables vanishes. 

Separating the mean and fluctuating parts of the unknown vector fields, transform 
(23) into an equivalent equation 

M{{u;o}v, {ho}., ^o) = - (Vx X (((vo)). x O - ((ho)), x (V, x H)), 

(((ho)).-V.)V-(((vo)).-Vx)H, -(((vo)). • V.)e) . (25) 

((vo)). and ((ho)), are independent of fast variables, and slow variables are not involved 
in the definition of the operator Al; thus, by linearity, a solution to equations (25), 
(8b), (8c) and (9) for n = can be expressed as 

iWoU {vo}., {ho}., ^o) = Co + E ((^o))/c + S^((/io)).) (26) 

fc=i 

(interpreted as an equality of 10-dimensional vectors). S^" can be determined from 
the equations 

v.xsr = sr, v.-sr = o, (sn.^o; (27a) 

V.xS^'" = S^'^ V.-S^'^'^O, (S^'> = 0. (27b) 

Sj^'(x, t) — {S'j^^, Sj^", S,f^, S'l^^) are solutions to auxiliary problems of type I, satisfying 
the boundary conditions similar to (2) and (4): 

MiS, >S,',5,')=(^— ,— ,— J, (28a) 

Vx • Sr = Vx • Sf = 0, (28b) 
(SD. = 0, (S^''% = 0, (28c) 

together with (27a); 
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A.(S^-Sr,^r) = (-V.x||,-|^,o), (29a) 

Vx ■ S^'^ = Vx ■ S^" = 0, (29b) 
(S^-'^), = 0, (S^''^), = 0, (29c) 

together with (27b). 

The problems (28), (27a) and (29), (27b) are particular cases of the problem (12); 
consistency relations (12d)-(12g) and solvability conditions (15) and (18) are trivially 
satisfied for them ( (15) is satisfied due to the assumption that f2 and H are globally 
bounded). Relations (28b), (28c), (29b) and (29c) hold true at any time, provided 
they do at t = 0. Therefore, any smooth vector fields, globally bounded together with 
their derivatives and satisfying (28b), (28c), (29b), (29c) and the boundary conditions 
similar to (2) and (4), can serve as initial conditions for the problems (28) and (29), 
as long as the respective solutions remain globally bounded in time. For CHM states 
V, H, ©, linearly stable to small-scale perturbations, the boundedness for any initial 
data is demonstrated in appendix B. Existence of steady or time-periodic solutions for 
generic perturbed CHM states, periodic in horizontal directions, is shown in section 
3. 

Equations (28a) and (29a) are equivalent to 

respectively. The second of these equations is an eigenvalue problem for eigenfunc- 
tions from ker A4' with non-zero spatio-temporal means of horizontal magnetic com- 
ponents. Since (0, (Ci,C2,0),0) G kei A4'* for constant Ci and C2, this problem has 
a solution regardless of whether kerA^ is trivial, as assumed in sections 5-10, or not, 
provided A4' does not have a Jordan cell associated with the zero eigenvalue. In 
computations it may be natural to consider, instead of vorticity components of (28a) 
and (29a), equations involving their vector potentials: 

/^vfov,v f^v,h qv,9 QV,p\ ^ * nv^ch^v ch,h nhfi cih,p\ un. 

^ Wfc > '^fe > "-"fc > "-"fe ) — K'^k ^'^k > "-"fe > ) — 

where {'V:>^Sfl^)h = 0, and to solve numerically auxiliary problems of type I in the 
terms of S^", S/', 5^'^ and S*^'^. 

^o(x, t, X, T) = (^0 , Co; Co) satisfics the equations 

^eo,eo,Ceo) = 0, (30a) 

Vx-C^ = Vx-CS = 0, (30b) 

(C). = 0, {to)h = 0, (30c) 

VxXCS = C Vx-C^ = 0, {Oh = 0. (30d) 

Relations (27), (28b), (29b), (30b) and (30d) guarantee that conditions (8b), (8c) 
and (9) hold true for n = 1. Again, it is sufficient to demand that (30b) and (30c) are 
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satisfied ai t — 0. Initial conditions for (30) can be found from (26) at t — 0. (Prom 
(24a), ((ho))h|T=o = (ho)ft|t=o-) Initial conditions for must belong to the stable 
manifold of the perturbed CHM state V,H, O, i.e. must exponentially decay in 
time. A permissible change of initial conditions for in the considered class implies 
a change of initial conditions for ^q, but then the respective changes in S^' and 
exponentially decay in time. 

6. Order and £^ equations: a-effect in the leading order 

Averaging the vertical component of (Al) for n = 2 over fast spatial variables and 
taking into account (13a), (8a)-(8c), (9) forn = and 1, (10) for n — 0, the boundary 
conditions for V, H, Vj and hj, and solenoidahty of Vq and ho, find 

^^"^'^^ + Vx X (V X (Vx X vo) - VVx • {vo}^ 



dt 

-H X (Vx X ho) + HVx • {ho}h)h = 0. (31) 
Substitution of the flow and magnetic field (26) now yields 



Here it is denoted 

r = (V X (Vx X - VVx • - H X (Vx X ^o') + HVx ■ ^0)1.; 

Kr,k = (V X (e^ X (Sr + e,)) - H x (e„ x S f) - \{Snm + li{Sl'')m)H, (33a) 
at^, = (V X (e^ X Si') - H x (e^ x (Sj''^ + e^)) - V(Si'')m + H(Om).; (33b) 

«Ii = (ali)i - «i)2 - {ahh, al, = (a^Ji, (33c) 

«1,2 = («l2)2, «2,2 = (42)2 - «2)l - «l)2, (33d) 

where (/)„ is the m— th component of a three-dimensional vector f; a^^, are also 
expressed in the terms of a^^^ by the last four formulae, where the superscript v is 
changed to h. The differential operator in the right-hand side of (32) represents the 
AKA-effect (anisotropic kinematic a— effect) operator in the leading order. 
Applying (17), obtain from (32) 

{0J2)v = {Mv (34) 



, fe=l m=l 



This expression is well-defined, if and only if the means ((/o fc '^^)) f*^^ m,k = 1,2 
are. This is the condition of insignificance of the AKA-effect in the leading order. It 
is stronger than the solvability condition (15a) for the system (A1)-(A3) for n — 2: 

= 0. (35) 
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A similar procedure reveals a possible presence of magnetic a— effect in the CHM 
system. Averaging of the horizontal component of (A2) for n = 1 over fast spatial 
variables with the use of (13b) and substitution of the flow and magnetic component 
of (26) yields 



i:(«H«o)>. + a2(Ao))*)+«'' ■ (36) 



.k=l 



dt 

Here it is denoted 

«^ = (V X Sf + {Sr + e,) X H),, = (V X (S^''^ + e,) + S^''^ x H),; (37) 

The differential operator in the right-hand side of (36) represents the magnetic 
a— effect. (17) applied to (36) yields 



(hi). = ((hi)),+VxX ( E {{[yidt}ivo}k + { l^c^Ut] (ho)) 




.k=l 

(38) 

Magnetic a— effect is insignificant in the leading order, if the means ((/q ot^dt)) exist 
and as a result l^i)h is well-defined by (38). The condition of insignificance of the 
magnetic a— effect is stronger than the solvability condition (15b) for the problem 
(A1)-(A3) forn = 1: 

io^D) = ((«!:)) = 0. (39) 

If the perturbed CHM state is steady or periodic in time, the 8 scalar relations 
(35) are sufficient for insignificance of kinematic a— effect, and the 4 scalar relations 
(39) are sufficient for insignificance of the magnetic a— effect. Note that vanishing 
of vertical components of the kinematic a— tensor in the vorticity equation or of 
horizontal components of the magnetic a— tensor in the magnetic induction equation 
is not required. For this reason we are speaking about insignificance and not about 
the absence of the a— effect. 

Thejterms {{/o I'^c?^^ in (34) and ^Jo^'^dtJ in (38) are not problematic, because 
and exponentially decay and thus the averages ((/q i^dtj and ((/q i^dt} are well- 
defined. It is shown in appendix D that ^Jo^'^dt]^ and i^Jo^^dt^ also exponentially 
decay. 

The multiscale approach remains feasible even, if (35) or (39) do not hold true. It 
is well-known (see DubruUe and Prisch, 1991), that in this case another slow time scale 
is appropriate, T — et. Then the new terms, d{u)i)v/dT and d{ho)h/dT, emerging 
in the left-hand sides of (32) and (36), balance the a— effect terms. (32) and (36) 
then constitute a closed system of equations (together with (8a) and (10) for n = 1). 
The mean-field equations turn out to be linear first-order PDE's; solutions to such 
equations generically exhibit unbounded exponential growth. Therefore, from now 
on we focus on a potentially more interesting case of an insignificant a-effect, where 
possible growth of perturbations may saturate due to nonlinearity. 
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7. Symmetries, gueiranteeing insignificance of the a— effect 
in the leading order 

In this section we consider the symmetries, in the presence of which the a— effect is 
insignificant in the leading order. They are compatible with the equations of thermal 
hydromagnetic convection and boundary conditions considered in this paper. 

A CHM regime O, V, H, is called parity-invariant with the time shift T, if the 
fields V and H are parity-invariant with the time shift T, and ft and © are parity- 
antiinvariant. A three-dimensional vector field f is parity-invariant with the time 
shift T, if 

f(-x,t) = -f(x,t + T), 

and parity-antiinvariant, if 

f(-x,i) = f(x,t + f); 
a scalar field / is parity-invariant with the time shift T, if 

/(-x,i) = /(x,i + f), 

and parity-antiinvariant, if 

/(-x,i) = -/(x,i + f). 

(We have assumed here that the origin of the coordinate system is located at the 
centre of symmetry on the mid- plane of the hquid layer.) If the perturbed CHM 
regime is parity- invariant, then we call a set of fields {uj,\,h,6) symmetric, if v and 
h are parity- invariant, and uj and 6 are parity-antiinvariant; it is called antisymmetric, 
if V and h are parity-antiinvariant, and lc and 6 are parity- invariant. 

We call a CHM regime ft, V, H, symmetric about the vertical axis xi — X2 — 
with the time shift T, if all the fields V, H, ft and are symmetric. A three- 
dimensional vector field f is symmetric about the vertical axis with the time shift 
f, if 

fl{-Xi, -X2,X3,t) = -fl{xi,X2,X3,t + T) , 
f2{-Xl, -X2, Xs, t) = -f2{Xl, X2, X3, t + f), 
fsi-Xl, -X2, X3, t) = f3{xi,X2, X3, t + f), 

and antisymmetric, if 

fl{-Xl,-X2,X3,t) = fi{Xi,X2,X3,t + f), 

f2{-Xl, -X2, X3, t) = f2{Xl,X2, X3, t + f), 
fsi-Xu -X2: X3, t) = -f3{xi, X2, X3, t + f); 

a scalar field / is symmetric about the vertical axis with the time shift T, if 

f{-Xi, -X2, X3, t) = f{xi, X2, X3, t + f). 
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and antisymmetric, if 



f{-Xi, -X2, Xs, t) = -f{xi,X2, X3, t + T). 

(Without any loss of generality the origin of the coordinate system is assumed to 
reside at the vertical axis.) When this symmetry is concerned, we call a set of fields 
(a;, V, h, 9) symmetric, if all the four fields are, and antisymmetric, if they are all 
antisymmetric. 

Symmetries without a time shift (for T = 0) are spatial, they can be possessed 
by CHM regimes of arbitrary time dependence. We will call such symmetries parity 
invariance and symmetry about the vertical axis, respectively. By contrast, only 
regimes, periodic in time with the period 2T can have spatio-temporal symmetries 
with a time shift T ^ 0. Travelling waves, for instance, can have such symmetries. 

Suppose the symmetry is defined by means of the operator S, i.e. symmetric and 
antisymmetric fields satisfy the conditions iSf = f and St = — f, respectively. Since 
all the symmetries considered in this section are of the second order, an arbitrary field 
f can be decomposed in a sum of a symmetric field (f + St)/2 and antisymmetric 
(f-<Sf)/2. 

If a CHM regime is symmetric, then symmetric and antisymmetric sets of fields 
are invariant subspaces for the linearisation operators C and Ai. In this case the 
right-hand sides of equations (28a) and (29a) are antisymmetric sets, and hence 
essentially S^' are antisymmetric sets (they are antisymmetric if initial conditions 
are; by construction, for any permissible initial conditions a symmetric part of a 
solution to auxiliary problems of type I exponentially decays). Consequently, (33) 
and (37) imply that in the presence of a symmetry the AKA- and magnetic a— effects 
are insignificant. Moreover, since 

«m,fe(^ + r) = -a^^kii) 
for a symmetric perturbed CHM state, the function 

^■{t)^{{l^a-^,,dt} 
satisfies a similar relation ■i?'(t + T) = — ■i? (t): 

rt i'2T rt ]^ p2T rt 

m = I '^m,k dt-^l I dt' di J. dt' di, 

whereby 

d\t) + ^-{t + T)^^j^ ^ a^^,{t') dt' + y. <,,(t') dt'^ dt 

lo + »m,k(t') dt'yi ^ 0. 
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8. Order equations: ct-efFect, insignificant in the leading order 

In this section we solve equations (Al)-(A3) ior n — 1, assuming from now on that 
the AKA- and magnetic «— effects are insignificant. 

The solution to equations (10) for n — 2, (8a) for n = 1 and (34) is 

-VxH + r, (40a) 

S Jo ^ ~ "^'^^ dx^dx, + ii L ^"-^ " "-^^ J ax^ax, J ' 

(40b) 

where the operator is inverse to the Laplacian in slow variables: for f (X), whose 
average over X is zero, g(X) = Vx^f is the mean-free solution to the equation 
Vxg = f , which is globally bounded together with the derivatives. In (40a) 

e = {{ /J ^'dt }} , r (X, t, T) = r - VxVi2(Vx ■ r). (41) 

exponentially decays in fast time (inheriting this property from ^q). It is shown in 
appendix D that also exponentially decays in fast time. 

Thus (vi)/i is expressed in the terms of ((fo))fc and ((/io))fc- Since the solvability 
condition (10) for the equation (9) for n = 1 is satisfied (see section 5), {^i}h can 
be found from (9). As shown in the previous section, vanishing of the average of 
(vi)/i over slow spatial variables guarantees that the condition (6g) for Vi|x=£(a;i,i;2) 
is asymptotically satisfied to any power of e. 

Transform the equations (Al)-(A3) for n = 1 using in (Al) the identity uJi — 
{uJijy + Vx X ((vo))fe: 

C'^i{u^l}v, {Vl},, ^hi}};,, 9i) + (Vx X i^0))h)3^^ - ((Vl), " V^)n 

+(((hi))^ • Vx)(V, X H) + 2i/(V, • Vx){u;o}. - V, x (H x (Vx x ho)) 
+Vx X (V X Wo + vo X 12 - H X (Vx x hg) - ho x (Vx x H)) 

+Vx x (vq X - ho X (Vx X ho)) + /3Vx6'o x 83 = 0; 

C\{Y^}h, {{hi}},) - ((vi);, ■ Vx)H + (((hiK ■ Vx)V + 2r/(Vx • Vx){ho}/. 

+Vx X (vo X H + V X ho) + Vx X (vo x ho) = 0; 

>C^({vi}ft, ^1) - ((vi), • Vx)e + 2/^(Vx • Vx)^o - (V • Vx)^o - (vo • Vx)^o = 0. 

In view of relations (26) and (40), (8) and (9) for n = 1, and due to linearity of 
this problem, it has solutions of the following structure: 
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fe=i 



dXm "^'^ dXra "^'^ dX^dX^dX, "^'^ dX^dX^dX, 

+ QZMivohivojm + Q^.^0)).((/l0))n. + Qt%{{ho)) dho} ) ■ (42) 

Here G;;; ^ = (G;;,";^, G^^;t, G^'^fe, J solve auxiliary problems of type II: 
mGZ>^) - -^m,k,sg^^ - 2u-^ - X (V X 
+(Sr + e^) X - H X (V. X Sf ) - Sf x (V. x H) + pSl'^e^) 



+Vx x (H x (e„ X Sr )) + /„ cxl,kdt j ^ (43a) 



i^m,k,j denotes the standard unit antisymmetric tensor), 
^"{GZk) = - e„. X (V X Sf + (Sr + e,) x h) + {{ aZ,udt }} 



(43b) 



£^(G-,) = -2«^ + + {{ jl <,,dt }} ^, (43c) 

V. X GZk = GT,k - em X sr, (43d) 
• Gj^'^fc = -(S'^''^)„i, Vx • Gj^y^. = ~('S'fe''')m, Vx • Gj^^fc = -(S'^''')„i; (43e) 

>C'^(G!:;;.) = -2^^ - X (V X S^''^ 
+S^'^ X 12 - H X (Vx X S^''') - (S^''^ + e,) x ( Vx x H) + /J^^^ea) 

>C^(Gi;,) = -2r;^ - x (v x (Sj'^ + e,) + x h) + {{ <,dt }} 



+Vx X (H X (e^ X (S^''^ + e,)) + f / a'^^.dt } —, (44a) 



m.i 

'^k 

/ n h h \ 11 ------\ ({ f U 



(44b) 

Vx X g;;;;, = g:;:;^ - e„, x s^^ (44d) 

Vx • gC;;^, = -(5,^"^)^, Vx • gC;:, = -(5^X, Vx • G^;;^, = -(5^'^)^; (44e) 
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Qm.fe = {Q'm% Q'm% Qm,fc> Q'm,k) so^ve auxiliary problems of type III: 

M-{QZ:k. QZ'l QZi) - Pm,kV^x ( - (sr + e.) X sr 

+Sf X (Vx X S;;;'') - (Sr + e^) X + S;;;^ X (Vx X S^,''')) , (45a) 

M\QZ:k, Q^t) = -Pm,.Vx X ((Sr + e.) X S;;;'^ + (Sr + e^) x Sf ) , (45b) 

M'iQZ'l QZi) = Pm,k (((Sr + e.) • Vx)^;;;^ + ((Sr + e^) • V^)Sl'') , (45c) 

Vx X Q-;,^ = Q™;- (45d) 



A^'^(Q;t' Qt', ) = -Vxx ( (sr + e.) x s::;'^ + s;:;'' x sr 

-Sf X (Vx X S;^'^) - (S;;;'^ + e^) x (Vx x Sf )) , (46a) 

M\Qt;k, Q^'i') = -Vx X ((S^,'" + e,) X (S:^^ + ej + S^;;'^ x Sf ) , (46b) 

A^'(Q;'r , Qtk) = ((Sr + e.) • Vx)^^^ + (S^^;^ • V^)Sl'', (46c) 

Vx X Q;^; = Q;^; (46d) 



-M^(Qrr,Q::;^c;fc) = p«.;^Vx x (-sr x s:r 

+{Si'' + e,) X (Vx X S:;;'^) - S^^^ x S^^ + {S^^ + x (Vx x Sj''^)) , (47a) 

M^Qt':,-, Qt't) = -Pm,kV^ X (S^ X (S;:;'^ + e^) + S;;;^ x (S^ + e,)) , (47b) 

M%Qt%-, Qt',f ) = Prn,k ((S^ ■ Vx)^^'' + (S^ ■ V^)S^/) , (47c) 

Vx X Q:;';;,'' = Q::,';^; (47d) 

Vx • Q^;^ - (48) 

{Pm,k = 1 for m < /c, pm,k = 1/2 for m ^ k, and Pm,k = =^ Q^'^ = Q^^^'^ = for 
m> k); Y'^ ,^ = 0^m,k,^m,k,^m,k,'>^m,k) solvc auxiliary problems of type IV: 

•'^H ^ m,ki ^ m,ki ^m,k) — Pm,k 







+ { J\al, - al,)dt }} ^ j (O, H, 9), (49a) 
Vx X y::;:, = y;;;- ; (49b) 
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_d_ 

dxi 



+ 



h,U! . 



(50a) 
(50b) 



Vx • Y- , = (51) 
(Y;; = for m > A;); and = (^^, {f) solve the problem 

C^iti, li, ti, Ci) = -2^(Vx • Vx)^^; + V. X (H X (Vx X ^o')) - Vx x (v x 

+eoXn-Uxi\/^x ^'i) - x (Vx x H)) - Vx x (vo x to + Co x i^o - C) 

-ho X (Vx X to) - ^0 X (Vx X (ho - ^H))) - /JVxCo X es, (52a) 
>C'(£L C?) = -277(Vx • Vx)£? - Vx X (£^ X H + V X £o') 

-Vx X (vo X to + £o X (ho - to)) (52b) 
J^'iCi, ^i) = -2«(Vx • Vx)eo' + (V • Vx)eo' + (vo • Vx)eo' + (Co ' Vx)(^o - ^o), (52c) 

VxX|^-|r = -Vxx|^, (52d) 

Vx-Cr = -Vx-C^, Vx-C^ = -Vx-C Vx-C? = -Vx-C (52e) 

Qmjfc' Y^ j^. and £i must satisfy the boundary conditions similar to (2) and 
(4). Means over fast spatial variables of their horizontal flow components and of 
vertical vorticity components must vanish, as well as spatio-temporal means of hori- 
zontal magnetic components. Clearly, the auxiliary problems are particular cases of 
the problem (12). 

Equations (45a), (46a) and (47a) are equivalent to 

fv({^vv,v f^vv,h (yjvfi ^vv,p\ _ ( (ov,v , _ \ ^ ov,u> 

+Sf X (Vx X S::;'^) - (Sr + e^) X Sr + S:;;'^ X (Vx X Sf )) , (53) 

^ \^m,k^^m,k^^m,k^^m,k) — ~\'^k + ^feJ X " X 

+Sf X (Vx X S:^'^) + (S:^'^ + e^) X (Vx X Sf ), (54) 

-St' X + (S:^'^ + e^) X (Vx X S^''^)) , (55) 

respectively, where (VxQm^fc)^ = 0- 

Together, (43e), (44e), (48), (51) and (52e) are equivalent to (8b), (8c) and (8e) 
for n = 1. It suffices to impose these conditions for the vorticity and magnetic parts of 
the solutions at t = 0, since the compatibility conditions (12e), (12f) evidently hold 
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true (for the auxiliary problems of type II this can be deduced from (28a), (28b), 
(29a) and (29b) ). Similarly, equations (43d), (44d), (45d), (46d), (47d), (49b), (50b) 
and (52d) are equivalent to equation (11) for n = 1. 

The spatial means of the initial conditions for the problems (43)-(52) can be 
determined averaging (16) over fast time: 

rt 



t=0 



-{{ e 



X 



Jo 



dt 



t=0 



-{{ e 



X 



Jo 



dt\, (56) 



t=0 W ^0 II h 

(the means exist by the assumption that the magnetic a— effect is insignificant); at 
any i > 

{^m,k)v = {Q.S,k)v = 0^m,k)v = {iTiv = {Q.m,k)h = 0^m,k)h = 0. (57) 

An initial condition for mean vorticity is determined by (24b): 



T=0 



t=0 



Averaging the magnetic part of (42) over spatial variables at i = 0, find 



-EE 



(hi) 



{^''m,k)h 



T=0 



(hi). 



t=o 



t=0 



T=0 



+ {^m'!k)h 



*=0 dXr,. 



T=0 



Now initial conditions for the problem (52) can be determined from (42) at t = 0. 

The choice of initial data must ensure that all solutions to the auxiliary problems 
are globally bounded with their derivatives. Generically, such solutions do exist, if 
the perturbed CHM state V, H, is periodic in horizontal directions and steady or 
time-periodic (see section 4); solutions are globally bounded for any smooth initial 
conditions, if the perturbed CHM state V, H, is linearly stable to small-scale per- 
turbations (sec appendix C). Modification of initial conditions for G^ ;,, Q^ ^ and 
^m,fc other permissible ones implies the respective changes in initial conditions for 
These changes must belong to the stable manifold of the perturbed CHM state 
V, H, 0, so that the resultant changes in solutions to the auxiliary problems and in 

decay exponentially in fast time. 

Since decay exponentially in time together with derivatives (see section 5), 
the right-hand sides of equations (52a)-(52e) also decay exponentially. As shown 
in appendices D and B for mean parts of the relevant components of ^^^d for the 
complementary fluctuating part, respectively, this implies an exponential decay of 
if the CHM state V, H, is linearly stable to small-scale perturbations. Similarly, 
under this condition, any changes in G^^^., Q^^^ and Y^^^ due to a permissible 
variation of the initial data for S^' exponentially decay. Otherwise, an exponential 
decay of must be ensured by an appropriate choice of the initial data. 

If the perturbed CHM regime has a symmetry considered in section 7, S^' are 
supposed to be antisymmetric sets of flelds and therefore the right-hand sides of 
equations defining auxiliary problems of types II-V are symmetric sets. Due to 
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invariance for the operator C of the subspaces of symmetric and antisymmetric sets 
in the presence of a symmetry, essentially G^j^., Q^ ^ and Y^ji. are symmetric sets 
(by construction the antisymmetric part of any permissible solution to the problems 
(43)-(50) exponentially decays). Then conditions (56) and (57) are automatically 
satisfied (except for vanishing of the mean of the vertical component of vorticity is 
not implied by parity invariance of the perturbed state). If the perturbed CHM state 
is steady or its symmetry is without a time shift, then the right-hand sides of auxiliary 
problems of type IV are zero and hence k~^- 

If the perturbed CHM regime is steady, periodic or quasiperiodic in fast time 
and/or spatial variables, it is natural to demand (in particular, for convenience of 
computation of spatio-temporal averages), that Sj^', G^ ^^, and j. are steady, 
or have the same periodicity or the same set of basic frequencies, respectively. Then 
and are decaying transients, bringing perturbations of the CHM system (in fast 
time) to their regular regime. 

9. Order and equations 

Averaging the vertical component of (Al) for n = 3 over fast spatial variables taking 
into account (13a), (8a)-(8c) and (9) for n = 0,1,2, the boundary conditions for 
V, H, Vi and hj and solenoidality of Vq and ho, find 

dt 9^ + ^Vx(^i). 

+ Vx X (V X (Vx X vi) - VVx • {vi}h - H X (Vx X hi) + HVx • {hi}^ 

+ Vo X (Vx X vo) - VoVx • {vo}h - ho X (Vx x ho) + hoVx • {ho}^)^ = 0. (58) 

Averaging this equation over fast time, substituting the flow and magnetic compo- 
nents of (26) and (42), employing (40a), (10) for n = 1 and (35), and recaUing that 
£0 and decay exponentially, find 

Vx X [-^((voK + z/V^((vo)). - (((vo)). • Vx)((vo)). + (((hoK • Vx)((ho)). 



2220/0/ 02 

+ dX.,((/^o))ijJ+ A:-,((t;o)).((t;ok+ Ar^^^^^^^ = 0, 

(59) 



10 1 

+ E E E ^ ^''r^,kA''^)t^ + D!:;l,((/io))fc + o^ o^ (f^^kA^^I^ 

ih.v 111 w \ \ I A vv,v II \\ II w I A 'vh,v II W // r, W i A hh,v 



where 



^llj = ((-^.g::;^^ - y{Gi:k)j + h,gi% + ii{Gi%),u (eoa) 

d:^:;, = ((-^.-^3, - y{YZh + H,Yl\ + H(y:'i),))„ (60b) 
GZk = G;:;^, + {{ [ al,^,dt }} e„ Yl^^, = Y^^^, - {{ [{a^^, - a^^,) dt }} e,; (60c) 







changing in (60) every occurrence of the first superscript "v" to "/i" , obtain expres- 
sions for D^^;^^^. and d^^;^^^.; 
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(61a) 

^m,k,j — \\ ^j^m,k * Wm,k)j ~r ^j^m,k "r "-Wm,k )j 

-iSn^St' - ('5^''').Sr + {Sl^^S^J^ + {S'J)jSI'))h, (61b) 



(61c) 

Averaging the horizontal component of (A2) for n = 2 over fast variables, substi- 
tuting (40) and flow and magnetic components of (26) and (42), recalling that and 
^'^ decay exponentially and taking into account relations (14), (39) and the boundary 
conditions for V, H, v„ and h„, find 

-^((ho)). + ^V^((ho)). + Vxx (| ((vo)), X ((ho)), 
+ E E {Dl>o}k + T^'^kihojk + a^^Vx^ Wi + d'liAhohfj 

+ A-ZiiMkiMm + A:!:iivoMho))m + AlJ^f ((/io)).((M..)) = 0, (62) 

where 

= ((V X g::;5, - H X g::;^,))., a;^^, = (v x y^I^, - h x y;;;:,)),; (63) 

changing here every occurrence of the first superscript "v" to "/i" , obtain expressions 

r -r\h,h j jh,h 

for ^m,k and d^^;.; 

AZk = ((V X QZi - H X Qr:r + Sr X S:;;'^))^, (64a) 

<i = ((V X Qj','^ - H X q;^;,^ + Sr X S^^;'^ + S^^;^ x Sf ))„ (64b) 

<'k = ((V X - H X Q':;::^ + S^^ x S:^'^)),. (64c) 



D''' are coefficients of the second order operators representing the so-called aniso- 
tropic combined eddy diffusivity. Simpler expressions for the operators arc presented 
in appendix E. d '' are coefficients of pseudodifferential operators, formally also of 
the second order, which can be regarded as representing an unconventional non-local 
anisotropic combined eddy diffusivity. All d '' vanish, if the perturbed CHM state is 
steady or possesses a symmetry without a time shift of a type discussed in section 7. 
To the best of our knowledge, such a physical effect was not encountered before. A "'' 
are coefficients of quadratic terms representing the so-called combined eddy advection. 

Equations (59) and (62) are solvabihty conditions for the system (Al)-(A3) for 
n = 2. To solve the order equations, determine from (58) an expression for (cc^s)^ 
similar to (34), then from (8a) and (10) for n = 2 deduce an expression for (v2)/j 
similar to (40), and finally solve the remaining parts of equations (A1)-(A3) for 
n = 2 (the fluctuating {-j^ part of (Al), fluctuating ^-^h part of (A2), and (A3) ) 
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to find expressions for {c«;2}i;, {^2}h, V^2^h and 62 analogous to (42). It is possible 
to construct in a similar way solutions to the systems arising for higher orders of 
and thus to obtain subsequent terms in expansions (7). 

(59) and (62) constitute a closed system of equations for mean perturbations (more 
precisely, for the leading terms in their expansions in power series in the spatial scale 
ratio). The vertical component of (62) and the horizontal component of (59) vanish 
identically. One can "uncurl" (59), introducing the large-scale pressure. However, 
since ((vo))^ and ((ho))^ are solenoidal in slow variables, it is more natural to solve this 
system numerically introducing stream functions for the two-dimensional flow ((vo))/i 
and magnetic field ((ho))/j and considering vertical components of (59) and of the curl 
of (62). 

Note that the Coriolis force enters only in the statement of the auxiliary problems 
and not in the mean- field equations (59). Derivation of these equations remains 
unaffected, if additionally the fluid is rotating in slow time, and the total angular 
velocity is re^ + e^T. Then the mean-field equation (59) for the flow involves an 
additional term representing the Coriolis force in the fluid layer rotating with the 
angular velocity r. 

10. Mean-field equations with a— effect terms: 

l£irge-scale perturbations of CHM regimes neeir a Hopf bifurcation 

As discussed in section 6, construction of the system of mean-field equations requires 
that the kinematic and magnetic a— effects are insignificant in the leading order. 
However, in the considered problem this implies the absence of a— effect terms in 
(59) and (62). In this section we consider perturbations of CHM regimes constituting 
branches parametrised by a small parameter, and show that if this parameter is linked 
with the scale ratio £, then the a— effect can remain insignificant in the leading order, 
but emerge in the mean-field equations. 

Let a family of CHM regimes admit an expansion in asymptotic power series: 

00 

(n, V, H, e) = ^ (f2„(x, t), V„(x, t), H„(x, t), 0„(x, t)) e\ (65) 

n=0 

^^n = Vx X V„. 

It can stem from similar expansions for the source terms F, J and S in (1). Perhaps, a 
more important example is a branch emerging from the CHM regime Vq, Hq, Go (not 
necessarily steady) in a Hopf bifurcation, occurring^ at /3 = /3o when two complex 
conjugate eigenvalues of L cross the real axis. Then the family of CHM regimes 
admits^ the expansion (65) for 

/3 = /3o + /32£' (66) 

^The quantity /3, proportional to the Rayleigh number, is chosen here as the bifurcation para- 
meter, because investigation of bifurcations in convection happening when this number is increased 
is a popular research subject (see, e.g., Podvigina 2006). Variation of any other quantity resulting in 
a Hopf bifurcation, e.g., r proportional to the Taylor number, can be considered; then the mean-field 
equations are not altered and there are minor modifications in the statements of auxiliary problems. 

^Unless the case is not generic, e.g., if the CHM state Vq, Hq, Oq has a large group of symmetries. 
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(see, e.g., Guckenheimer and Holmes 1990). The deviation of the emerging CHM 
regimes from the one at the critical point 13 — (3q \s asymptotically close to a sum 
of eigenvectors associated with the pair of imaginary eigenvalues of M. at this point, 
i.e. r2i,Vi,Hi,0i belongs to the subspace spanned by these eigenvectors (which is 
generically two-dimensional) . 

We assume (66), where positive and negative correspond to a direct and reverse 
Hopf bifurcation, respectively, or, for j32 = 0, the expansion (65) is due to dependence 
of the source terms F, J and S on e. Everywhere in this and next section we refer to 
the operators £ and A4 in which fig, Vq, Hq, Qq and Pq are substituted in place of 
f2, V, H, 9 and /? in (5) and (20). The same substitution is assumed in any equation 
and definition of any quantity introduced in previous sections and referred to here 
or in the next section. Large-scale perturbations satisfy (6) and are sought in the 
form of the series (7). Construction of the mean-field equations differs from that for 
£— independent CHM regimes only in a slight complication of algebra: new terms 
involving higher-order terms fin, V„, H„ and 0„ for n > 1 emerge in the right-hand 
sides of the equations to be solved. 

Unless specifically mentioned, the CHM regime Vo,Ho,©o the point of bifur- 
cation is not required to be symmetric. If it has a symmetry considered in section 7, 
then symmetric and antisymmetric sets of vector fields constitute invariant subspaces 
of A4, and generically Qi, Vi, Hi, 6i is either a symmetric, or antisymmetric set. 

Substitution of (7) into (6) gives rise to equations (A4)-(A6) (see appendix A). 
Using (A4) for n = and (A5) for n = 1, (24) can be verified. The equations at order 
£° coincide with (23). Their solutions have the structure (26), where vector fields 
S^'(x, t) satisfy (27) and solve the auxiliary problems (28) and (29), and exponentially 
decaying mean-free transients ^o(^' ^) solve (30). If the CHM state Vq, Hq, Qq 
possesses a symmetry of a type considered in section 7, then S'^' are antisymmetric 
sets, and a-effect is automatically insignificant in the leading order. 

Averaging of the vertical component of (A4) for n = 2 yields (31) and, conse- 
quently, (vi)/j is determined by (40), provided the a— effect is insignificant in the 
leading order. Therefore, solutions to equations (A4)-(A6) for n = 1 have the struc- 
ture 



({u;i}„, {vi},, ih,u ^i) = + E SkMk + sl'-iMk + si'-ivoh + sl'-{{hoh 



I'QIlm 




m=l 



+0'^^'-((h)) ((h)) I v^,- ^'Vx , y/.,. g^^Vx ((Ml ^ 



(67) 



Here j^., Q^';;. and are solutions to the auxiliary problems of types II-IV; 

S;^' = (S^'^, S;^", S;^'', S;^^) satisfy the boundary conditions of the kind of (2), (4) and 
solve auxiliary problems of type V: 
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= (- Vx X (Vi X Sr + (Sr + e.) X - Hi X (V. x Sf ) - Sf x (V. x Hi)), 
-V. x ((Sr +e,) x Hi + Vi x Sf ), (Vi • Vx)5,^'^ + ((S^ + e,) • Vx)ei) , (68a) 

Vx X sr = sr, Vx ■ sr = Vx ■ sr = Vx ■ sr = o; (esb) 

(SD. = {^I'lh = {Sl'lk = 0, (68c) 

Ai(sr,sr,^r) 

= (- Vx X (Vi X sr + sr X Oi - Hi X (Vx X sr) - (Sj''^ + Bk) X (Vx X Hi)), 

-Vx X (Sr X Hi + Vi X (Sr + e^)), (Vi • Vx)^?'' + (Sr • Vx)ei) , (69a) 

Vx X sr = S^'^ Vx ■ S^'" = Vx ■ sr = Vx ■ sr = O, (69b) 
(Sr). = = (Sr)/. = 0. (69c) 

Solvability (compatibility) conditions (12d)~(12f) and (18) for the problems (68) and 
(69) are easily verified. 

Averaging of the vertical component of the vorticity equation (A4) for n = 3 over 
fast spatial variables yields 



dt dT +^^x(^i). (70) 

+Vx X (Vo X (Vx X Vi) - VoVx • {Vi};, + (Vi + Vo) X (Vx X Vo) - (Vi + Vo) Vx • {vo}fe 

-Ho X (Vx X hi) + HoVx ■ {hi},, - (Hi + ho) x (Vx x ho) + (Hi + ho) Vx ■ {ho},,),, = 0. 

Averaging (70) over fast time upon substitution of (67), obtain an equation for 
the mean flow perturbation, generalising (59): 

Vx X (-^((vo)). + ^V^((vo)). - (((vo)). • Vx)((vo)). + (((ho)). • Vx)((ho)). 



h,v 



j=l m 

^ ^ 1 



^ dXudX^ (d^,fej((^o))i + dm,fej\\"'0//i 



+ K::hi^oMvo))m + K::UvoMho))m + A:jr((^0))r(/^0))n.jjj = 0. (71) 

Here D^^^, d^^^j and A^'^'^j are determined by (60) and (61), 

A^:, = i-ivo), sr - Vo (^r ), - (vo, (sr + e„) - vi (^r).- 

+(^o).- S;^'^ + Ho {Sf)m + Sl^ + Hi (5^'^),))., (72a) 
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+(i/o), S::;'^ + Ho {S^^'U + {Stt + em) + Hi {S^'),}n. (72b) 

Averaging (A5) for n = 2 over fast variables upon substitution of (67), obtain an 
equation for the mean magnetic field perturbation, generalising (62): 

-^((ho)). + ^V^((ho)). + Vx X (((vo)). X ((ho)). + E(^f ((V)fc + ^^'((Mfc 

+ K:iivoMv,)U + <iivo))kiho)U + Al^if ((/.o)).((/^o))n^))) = 0. (73) 
Here D^^^ d^ ^^^^- and A;;;;^ are determined by (63) and (64), 

a;/' = X Ho + (Sr + e,) X Hi + Vo X Sf + Vi x Sf ))., (74a) 
Al^"" = ((S^'^ X Ho + S^'^ X Hi + Vo X S^''^ + Vi X (S^''^ + e,))),. (74b) 

The new terms involving coefficients .4^j and >A{^ represent the AKA- and 
magnetic a— effects, respectively. If the CHM regime Vo, Hq, 0o is symmetric, but 
f2i, Vi, Hi, ©1 do not constitute a symmetric set (e.g., if a symmetry-breaking Hopf 
bifurcation is examined), then S'-' have non- vanishing symmetric parts and generically 
the coefficients (72) and (74) are non-zero. Thus, like in the almost axisymmetric 
models considered by Braginsky (1964a-d, 1967, 1975) and Soward (1972, 1974), 
a small deviation of the perturbed CHM regime from symmetry gives rise to the 
a— effect in the mean-field equations. The combined eddy diffusivity operator in (71) 
and (73) can be simplified using the identities of appendix E. 

11. Mean-field equations with the a— effect and cubic nonlineeirity: 
Icirge-scale perturbations of CHM regimes neeir a pitchfork bifurcation 

In this section we consider again a family of CHM regimes represented by the power 
series (65) in a small parameter, which we link with the ratio e of spatial scales 
in large-scale perturbations like in the previous section. We assume now that the 
operator M. (considered in this section, as C and any equation or quantity introduced 
previously, for f2 = f2o, V = Vq, H = Ho, © = ©o and P — Pq) has a non-trivial 
kernel including a small-scale eigenvector S'(x, t) — (S'^, S^, S'*, S^): 

M{S'^,S^,S^) = 0, VxXS^ = S'', 

Vx ■ = Vx ■ = Vx ■ S" = 0, (S-), = (S'^), = (S'O^. = 0. 

Generically, when exists, S' spans keiAi. (If the CHM regime Vo,Ho,0o is non- 
steady, an eigenvector in the kernel of is a globally bounded in space and time 
solution to the system of equations (19), (6d)-(6g), which does not decay in time.) 
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Like in section 10, we assume (66) and investigate stability of CHM regimes, whose 
dependence on e (65) is due to a similar dependence on e of the source terms F, J 
and S in (1) (for = 0), or a family of CHM regimes emerging from Vq, Hq, Oq (not 
necessarily a steady state or a periodic orbit) in a pitchfork bifurcation at /3 = /3o, 
occurring when the Rayleigh number is varied and an eigenvalue of Ai passes through 
zero. A positive or negative (52 corresponds to a direct or reverse bifurcation, respec- 
tively. The family of CHM regimes emerging in the pitchfork bifurcation admits^ the 
expansion (65) (see, e.g., Cuckcnheimer and Holmes 1990). Their deviation from the 
regime at the critical point /3 = /3o is asymptotically close to the mean-free eigenvector 
from ker Al: 

(12i, Hi, Vi, Oi) = x(S^ S", S^ S^). (75) 

Construction of mean-field equations differs from the previously considered cases 
in the following aspects. Now, a new solvability condition for the system (21) must 
be taken into account: orthogonality of the right-hand side of the equations (21a)- 
(21c), derived from (A4)-(A6) for n = 2, to the eigenvector S;^(x, t) G kerA^*. 
Also, expressions for the first terms in expansions of solutions consist now of more 
terms, involving solutions to a larger number of auxiliary problems. Insignificance 
of the a— effect in the leading order incorporates now more conditions than we have 
derived in section 6, and it is unlikely that they all can be satisfied by parameter 
tuning. Under these circumstances, in this section we impose symmetries of the 
leading terms in the expansion of the perturbed CHM regime (65). We assume that 
the CHM regime Vq, Hq, ©o the point of bifurcation has a symmetry considered in 
section 7. Consequently, symmetric and antisymmetric sets of vector fields constitute 
invariant subspaces of M.. We assume in addition that S (x, t), and hence S;^(x, t) 
and r2i,Vi,Hi,0 1 arc antisymmetric sets. Then the a— effect is insignificant in the 
leading order and all solvability conditions for auxiliary problems are automatically 
satisfied. 

Upon substitution of (7) into (6), the leading terms of the resultant equations, at 
order e'^, yield equations (23). Their solution can be expressed as 

(u;o, vo, ho, ^o) = €o + E {^k{Mk + Sj'-((/io)).) + co(X, T)S- + (0, ((vo)),, ((hoK, 0). 

k=l 

(76) 

Here S;^. (x, t) arc antisymmetric sets satisfying (27)-(29), and exponentially decaying 
mean-free transients ^q(x, t,X, T) solve (30). 

Higher order equations (A4)-(A6) in the hierarchy are the same as in the previous 
section. Averaging the vertical component of (A4) for n = 2 obtain (31); upon 
substitution of the fiow and magnetic components of (76), find from (10) for n — 2 
and (8a) for n = 1 

(vi). = ((v,)). + r - VxH + i: ^ {{ f <,dt }} ^ 

fc=lm=l \ "'O iJ^m 

''See footnotes 5 and 6. 
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where and a!^ f, are determined by (33), 

2 



2 /-"^ S^Vx' 
ti „ti \Jj.ll -X. 



m=l 



Equations (A4)-(A6) for n = 1 have a solution 

dco 



(a;i, vi, hi, ^i) = ^1 + S-ci(X, T) + Sco + Q'^'^'c^ + ( G^'' 

fc=i \ 



+ Qm:kiho}kiho}m + ^Zkgx,dX^dX, ^ ^ '^''^ dX^dXr^dX, 

Here G;^;^, 0^;;;^, Q;;";^, Q;;(';;:,, Q|;J^, Y;;;';,, Y^;^ and S;^' are solutions to the auxihary 
problems of types II-V; G^'' = (G^''^, G^'^, G^''^, G^'^) solve auxiliary problems of type 

£-(G^'-) = -2z/- efe X (Vo X S'^ + S'' x 12o - Hq x (V^ x S'*) 



2 



-S'^x(VxXHo)+/3o^''e3)-VxX(Hox(efcXS''))+E{/ (^9^) 



m=l -0 ' ^ dx^ 

£'^(G^'-) = -2r;^ - e, X (Vo X S'^ + X Ho) + E { / }} (79b) 

= -2.|^ + (yo).^^ + E {{ r }} (79c) 

V.xGr = G^-efeXS^ (79d) 

V. • Gr = -(^")fc, V.-Gr = -(5% V. • Gf = -(5^^)^; (79e) 
Q^v = (Q^ ''^, Q^-'", Q^-''^, Ql^) solve auxiliary problems of type III': 
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M{Qr,QT',Qr) 

= (Vx X (-(S^'" + Bk) X S'^ - S'' X Sr + Sf X (V, X S^) + S'^ X (V^ x )), 

-V. x ((Sr + e^) X S'^ + X Sf ) , ((Sr + e,) • V^)S' + (S^ • V^)Sl'') , (80a) 

V. X Qr = Qr ; (80b) 

A^(Qf'",Qf'',Qf'') 

= (Vx X (-S^''' X S'^ - X S^'" + (S^'^ + e,) x (Vx x S'') + x (Vx x S^'^)), 

- Vx X (S^'^ X S'^ + S'^ X (S^''^ + efe)) , (S^'^ • Vx)^' + (S'' • Vx)^^''') , (81a) 

Vx X Qf = Qf'-; (81b) 

MiQ^^'"", Q'"''', g"'^'^) = (Vx X (-S^ X S'" + S'* X (Vx x S'^)), 

-Vx X (S^ X S'^), (S'' ■ Vx)^^) ; (82a) 
Vx X Q'='='^ = Q-^^'^; (82b) 

Vx ■ = Vx • Qr^' = Vx ■ Qr' = 0; (83) 
^rn,k,j solve auxiliary problems of type IV': 

A<(Y^^,^, Y^;^,^, r^l^.) = - {{ }} -^{Qo, Ho, eo), (84a) 

Vx X y::^,^^ = y;;^,^, Vx • Y;;;;,,^. = O; (84b) 
S' = (S'^, S^, S'^, S^) solve auxiliary problems of type V: 

A4(S^,S^5^) = (Vxx(-VixS^-S"xrii +Hi X (Vx x S'*) + S'^ x (Vx x Hi)), 

-Vx X (S'' X Hi + Vi X S'^), (Vi • Vx)^' + (S^ • Vx)ei) , (85a) 

Vx X S'" = S^ Vx • S'" = Vx • = Vx • S'^ = 0. (85b) 

Solutions to these auxiliary problems must satisfy the boundary conditions of the 
kind of (2) and (4), and have zero means over fast spatial variables of horizontal 
components of flow parts and of vertical components of vorticity parts must vanish, 
and spatio-temporal means over fast variables of horizontal components of magnetic 
parts. Due to the assumed symmetry of the CHM regime (Vq, Hq, ©o) and antisym- 
metry of the sets of fields (Vi, Hi, 0i), S' and S^', the right-hand sides of the systems 
of equations (79)-(85) are symmetric sets. Thus all solvabihty conditions for these 
problems are satisfied, and the solutions Q*^ '', ^m,k,j symmetric sets. 

Averaging of the vertical component of the vorticity equation (A4) for n = 3 over 
fast spatial variables yields (70). Averaging it over fast time and substituting (78), 
obtain an equation for the mean flow perturbation, analogous to (71): 
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Vx X (^-^((voK + ^V^((voK - (((voK ■ Vx)((voK + (((ho)), ■ Vx)((ho)). 
+ E ^ f-^r^o + + E ( ^ ( E ( D- ,^.((-o)). + ^'r^lAhoh 

j=l \ m=l \ ^^rn \ \ 

-Vx (dm!fcj((^^o))l + d^lj{{ho}l) + Edm!fej,i QX^X, I + ^"r^J^O ] 
v,v II \\ I Ah.v II J \\ , A cv,v II w , A ch,v , 



+'^'',jivo}m + ^j((/io))m + AZjivojniCo + A^J^'] {ho} ^Cq (86) 

+ E (a;X((^o)).((^o))^ + AiV((^o)).((/io))^ + At%iho))Aho))rn)^^ ) = 0. 
Here DX,fc,i) d^^^^j) A^^^j and are determined by (60), (61) and (72), and 

^Z, = ((-(^o),G^^ - Vo(G^''), + (i/o),G:;^ + UoiC^^)),. (87a) 

= ((-(H)^YX,,, - Vo(y^':,,). + (^o).y:;^, + Ho(y^l,).K, (srb) 

G;;'' = G;;^ + E { [ <,k dt }} e„ Y^;^,^^ = Y^^^,^^ - {{ f a^,^, dt }} e,; (87c) 



k=i 

= ((-(^o).-S^ - Vo(^^), + {Ho),S^ + Ho(^''), 

-{V,),S^ - V.iS^j + {H,)jS'^ + Il,{S%))f,; (87d) 

A,^ = ((-(H)iQ'"'" - yoiQ^'^j + (^o)iQ'"'" + VLo{Q'''% - S^iS^)j + S\S%))^; 

(87e) 

AZ^ = i-{Vo)jQZ'' - y.iQZ'io + w^q:::"" + Ho(gj.^''^), 

-(5''),sr - s^(^r).- + (^').s;;;'^ + s'^(5^'^),)),, (87f) 

A^'J = ((-(^o).-Q^''" - ^oiQtnj + (^o).-Q^'''^ + Ho(g^^''^), 

-{S^),St' - S^(5^"), + (5'^),S::;'^ + S\S^^%))n. (87g) 

Averaging (A4) for n — 2 and substituting (67) obtain an equation for the mean 
magnetic field perturbation, analogous to (73): 

-^((ho))/. + r?V^((ho))/^ + Vxx ((vo)), x ((ho)), + A'^^ + A^ 
+ E f^f ((^o)). + ((/^o)). + Df ^ + Ar''^((^o)).co + Af''^ihohco'j 

+ E E (^l;; ( d;:;!^((^o)). + j^'j:Ahoh + e 

(d:;;!^((^o))i + d:;;;^,((Mi)) (§8) 



fe=i 



A;t^o)).((^o))n. + A;^^ ((^o)).((/^o))^ + Ai%\{hoMho}m 11=0. 
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Here 'D'^l^f.j, d';!^,, and A;;;^^ are determined by (63) and (64), A^" by (74), and 
Df = ((Vo X G^''^ + G^'^' X Ho)).; d:^,, = ((Vo x y;;;^,^. + Y^^^^,^^ x Ho)).; (89a) 

A"'^ = ((Vo x S'* + Vi x S'* + S'' X Ho + S" X Hi))„; (89b) 

A'^ = ((Vo X Q'^^''* + Q'^^'" X Ho + S'' X S'*))^; (89c) 

A^''''^ = ((Vo X Q™''^ + Q™'" X Ho + X + S^''' x S'^))„ (89d) 

Af -'^ = ((Vo X Qf ''^ + Qf'"" X Ho + X S^''* + S^'" x S'*)).. (89e) 

If a family of CHM regimes emerging in a pitchfork bifurcation is considered, com- 
parison of (82) and (85) with the use of (75) reveals that S' = 2xQ^'^''; consequently, 

= 2xA] and a'' = 2xA^ . 

Equations (86) and (88) constitute a closed system together with the equation 
for the amplitude Cq, which is the solvability condition for the equations (21a)-(21c) 
derived from (A4)-(A6) for n = 2, consisting of orthogonahty of their right-hand 
sides to the eigenvector ^\ e kcrTW*. As before, orthogonality is understood in the 
terms of the scalar product ((a' -b )), where a' = {a'^,a'*,a^} and b' = {b^,b'^,6^} 
are 7-dimensional vector fields in the layer. 

To derive this equation we need to calculate ^'ijh- Using (24b), find from (70) 

(wg). = ((wg)). + (Vxx {{( Vo X (Vx X Vi) - VoVx • {vi}^ + Vi X (Vx X Vo) 

-ViVx • {vo}/, - Ho X (Vx X hi) + HoVx ■ {hi},, - Hi x (Vx x ho) + Hi Vx ■ {ho}/, 

+Vo X (Vx X Vo) - voVx • {vo}/, - ho X (Vx x ho) + hoVx • {ho}^ di}} . 

Substituting expressions (76) and (78) and "uncurling" this equation, using (10) for 
n = 2, obtain 

(V2),. = ((V2)),. + # - Vx(Vi2(Vx • $)) + 6, (90a) 
where ^2 decays exponentially in fast time and 

*(X, T, t) = f b^'^c, + A]cl + E f ^ (d- CO + E f ,,.((^0)). 

j=l \ m=l \ \ fc=l \ 

+D::;:,,,((/io)). + ^..((^o))! + ~^t\A^^)i^) + e<^.,. Sfe)) 

+ Xlj((^^o))m + ^mj((/io))m + A^'J ((vo))mCo + At']{ho)lmCo 
+ E [^Z%A{voMv.}rn + ^tUv^^^^^^^ . (90b) 

The coefficients in (90b) are defined by the modified formulae (60), (61), (72) and 
(87) for the respective (without tildes) coefficients of the equation for the mean per- 
turbation of the flow (86), where the following change is implemented: A coeflicient 
in (86) is a spatio-temporal mean of the horizontal component of a certain vector 
field, say, C(x, t); then the respective coefficient in (90b) is equal to ^Jo^{C)h^ dt^. 
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If the symmetry of the CHM state Vq, Hq, ©o is without a time shift, then $ = 0; 
otherwise the symmetry properties of various vector fields assumed in this section 
imply that # is T— periodic in the fast time. 

Consider the operator Ai, the restriction of Ai' onto the subspace, defined by the 
condition ((h)),, = 0. The adjoini^operator is M* = {{M'*f, UM'*)^}h, {M'*^) (its 
domain is the same as that of M.). If S;^(x, t) e ker Al*, then 

S, = (s:^, S^- {{ / |((V X S:) X (V X H) -V X (V X S':))^} dt } , S^^ e kevM*. 

It is more convenient to employ S;^, than S'^. We normalise^ S;^ so that ((S;^ • S )) = 1. 
Substituting 

(u;2,V2,h2,^2) = {u:',y,X,02) + E(SrWfc + St((/l2))fc) 

k=l 

+ {'^^x{^l)h,{^2)h,i^2}h,0), (91) 

(76)-(78) and (90) into (A4)-(A6) for n = 2 and scalar multiplying the result by 
S'^, obtain the equation under discussion. It docs not involve any unknown quantity 
except the mean fields ((vo))/i and ((ho))/j and the amplitude Cq for the following reasons: 
((£(u;'2, v^,h'2,^^) ■ S,)) = 0, since (w'^), = (v^), = ((h^)), = and S,(x,t) e kerM*; 
(('f2))fe and ((/i2))fc do not appear, since vector fields Sj:'(x, i) are solutions to auxiliary 
problems of type I; all terms involving ((fi))fc, ((^i))fc, Ci and their derivatives do 
not contribute, since by our assumptions r2o,Vo,Ho, Oo is a symmetric set, and 
r2i,Vi,Hi,0i, S' and S;, are antisymmetric ones, and hence these terms constitute 
symmetric sets {{vo}^, iho}k, Cq and their derivatives enter into the expressions 
for u?0)Vo,ho,^o via terms, which are antisymmetric sets, and the expressions for 

vi, hi, ^1 via terms, which are symmetric sets; opposite to this, ivi}k, ((^i))fc and 
Ci enter into the expressions for a;i,vi,hi,^i via terms, which are antisymmetric 
sets); terms involving are absent because these vector fields and their derivatives 
exponentially decay in fast time. Thus, when deriving the equation one can set ^„ = 0, 
((vi)),. = ((hi));, = ((v2))^=((h2));, = and ci = in (76)-(78), (90a) and (91). 

The resultant equation is an evolutionary equation for cq. (It is bulky, and for 
this reason it is presented in appendix F.) If the CHM state Vq, Hq, Qq does not have 
a symmetry without a time shift, linear and nonlinear pseudodifferential operators 
involving the inverse Laplacian are present in the equation - for instance, all operators 
employed in Vx(Vx^(Vx • ^))- The terms, quadratic in Vq, Vi, ho and hi in (A4)- 
(A6) for n = 2 give rise to cubic nonlinearity. 

^Suppose for simplicity that M.* is an elliptic operator not having Jordan cells of size 2 x 2 or 
more, A1*s*p= ApS*p. The eigenvectors s*p constitute a complete set. Let us expand S': 

oo oo 

S'=^a-pS*p =^ (|S'p) = ^ (Tp(s*p • S') ^ (s»p • S') 7^ for somep. 
p=i p=i 

However, scalar multiplication of the eigenvalue equation for s*p by S' yields = Xp{s^p ■ S ), since 
S- G ker A^. Thus (s*p • S ) only for s,p e keiM*, and hence (S;, • S') = (S;, • S') ^ 0, since 
(S'')^ = 0. 
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12. Concluding renicirks 

We have considered weakly nonlinear stability of small-scale Boussinesq CHM regimes 
to perturbations, involving large spatial and temporal scales. Equations for the mean 
parts of the leading terms in power series expansions of perturbations in the scale ratio 
have been derived, using homogenisation techniques. In general CHM regimes exhibit 
the a— effect in the leading order; consequently, the mean-field equations turn out to 
be linear first-order partial differential equations, which generically have exponentially 
growing sohitions. If the a— effect in the leading order is absent (e.g. if the perturbed 
CHM regime is symmetric about the vertical axis or parity-invariant with or without 
a time shift), then the resultant mean-field equations (62) and (59) are nonlinear 
second-order partial differential equations, generalising the standard equations of 
magnetohydrodynamics. They possess new terms, describing eddy diffusivity and 
eddy advection, both anisotropic. If the CHM regime is non-steady and lacks a— effect 
in the leading order not due to a spatial symmetry without a time shift, new physical 
phenomena emerge, which are described by non-local pseudodifferential operators. If 
the regime is close to a symmetric one, e.g., it is near a symmetry breaking bifurcation, 
the a— effect emerges in the mean-field equations as well. Close to the point of a Hopf 
symmetry-breaking bifurcation the mean- field equations take the form (71) and (73). 
Near a pitchfork symmetry-breaking bifurcation, the mean-field equations (86) and 
(88) involve a new scalar amplitude of an additional mean small-scale neutral mode 
of the operator of linearisation. An equation for the evolution of this amplitude, 
presented in Appendix F, has a cubic nonlinearity. If, in addition, the perturbed CHM 
state is non-steady and the symmetry in the system is spatio-temporal (with a non- 
zero time shift T), new terms involving a non-local operator (the inverse Laplacian) 
acting on quadratic products of the mean perturbation components emerge in the 
amplitude equation. 

No boundary conditions have been set for the mean fields of perturbation in the 
horizontal directions, since they do not affect our derivation. They can be imposed at 
infinity, or on the boundary of a finite region in slow coordinates (e.g., of a region, the 
size of which grows as e^^ in the original fast coordinates); they must be chosen to 
serve specific applications of the theory that has been developed here. An important 
requirement is that the perturbation remains globally bounded - otherwise the basic 
assumption on the smallness of the perturbation and the asymptotic nature of the 
expansion become broken. 

In general, the mean-field equations do not satisfy the usual equations of energy 
balance, suggesting that their solutions are not uniformly bounded in time, and may, 
perhaps, collapse in a finite time. Such a growth, of course, makes the equations 
unsuitable for description of the subsequent evolution of the perturbation. For in- 
stance, the evolution of a perturbation of a CHM regime emerging in a pitchfork 
bifurcation is significantly affected by the sign of the coefficient in the term Ccq 
in the equation for the amplitude of the neutral mode: < damps the growth 
of the amplitude and stabilises the perturbation, and C"^ > supports its infinite 
superexponential growth. 
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We have derived the mean field equations for a perturbation of a CHM system 
in a horizontal plane layer, translation invariance in which is broken by the presence 
of source terms in the governing equations (l)-(3). Amplitude equations for a large- 
scale perturbation of a source-free CHM system will be derived in a sequel to the 
present paper. In the absence of the sources, or if an unsteady CHM state emerging 
under the action of steady sources is considered, the structure of the null space of the 
operator of linearisation Ai is significantly different from the generic one, assumed 
here, resulting in a different set of amplitude equations (not all of which will be 
evolutionary). In a future study we will also analyse numerically solutions to the 
amplitude equations for steady and periodic in time CHM regimes, symmetric about 
a vertical axis. In particular, this will reveal the nature of physical effects described 
by the new operators emerging in the amplitude equations. 

A similar derivation can be performed for a perturbation of a CHM regime in a 
spherical shell by considering the equations of hydromagnetic convection in spherical 
coordinates. These equations will be more appropriate for geophysical applications, 
but their structure is significantly more complex. 
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Appendix A. The hierarchies of equations 
for weakly nonlinear perturbations 

The following equations arise at order after substitution of the series (7) for the 
perturbation into (6a)-(6c) and expansion in power series, for the perturbed CHM 
state f2, V, H, © independent of e: 

V„, h„, On) - + U (2(Vx • Vx){Wn-l}. + V^W„_2) 

+Vx X (V X UJn-1 + Vn-1 X O - H X (Vx X h„_2 + Vx X h„_i) - h„_i X (Vx X H) 

n-2 

+ Xl(Vfc X UJn-2-k - hfe X (Vx X hn-2-k + Vx X h„_3_fc))) 
fe=0 

ri— 1 

+Vx X (^(Vfc X UJn-l-k - hfe X (Vx X h„_i_fe + Vx X K-2-k)) 
k=0 

-H X (Vx X h,_i)) + /?Vx^n-i X es = 0, (Al) 
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£'Xv„, h„) - + r/ (2(Vx • Vx){h„_i}, + V^h„_ 

n—2 n— 1 

+Vx X (v„_i X H + V X h„_i + ^ Vfc X K-2-k) + Vx X ^ Vjfc X K-i-k = 0, (A2) 

k=0 k=0 

^\^n, 9n) - ^ + ^ (2(Vx ■ Vx)^n-1 + V^^.-a) 

n— 1 n—2 

-(V • Vx)^n-1 - E(Vfc ■ Vx)^n-l-fe - EK ■ Vx)^n-2-fc = 0. (A3) 
fe=0 k=0 

The following equations result from (6a)-(6c) at order £" for the perturbed CHM 
state (65) for the parameter dependence (66): 

£-(u;„, Yr^, K, 9n) - + (2(Vx • Vx){a^n-l}. + 

n-1 



+Vx X (E C^fe X Wn-l-fe + V„_i_fe Xftk-HkX (Vx X X (Vx xHfe)) 

k=0 

n-2 

+ E(Vfe X <^n-2-k - Hfe X (Vx X h„_2-fe) - hfe X (Vx X h„_2-fe + Vx X h„_3_fe))) 
fe=0 

n 

+Vx X (EC^fe X U>n-k + ^n-k X fife - Hfe X (Vx X hn— fe) hjj_fe X (Vx X Hfe)) 



k=l 

n-1 



+ E(vfc X ujn-i-k - Hfe X (Vx X h„_i_fc) - hfc X (Vx x K-i-k + Vx x h„_2- 



fe=0 



+ (/92Vx^n-2 + /32Vx^n-3 + /^oVx^n-l) X 63 = 0, (A4) 



£'^(V„, h„) (2(Vx- Vx){hn_i}ft + V^h„_2) + Vx X E (Vn-fc X (Hfe + hfe_i) 



k=l 

n-1 



+Vfe X h„_fe) + Vx X E (vn-i-fc X (Hfe + hfc_i) + x h„_i_fc) = 0; (A5) 



k=0 



an n 
jC-'iVn, On) - + K (2(Vx • Vx)^n-1 + V^^n-s) - E ((^fe ' Vx)^n-fc 



fe=l 



+ (v„_fe ■ Vx)(efe + Ok-l) + (Vn-k ■ Vx)efc_i + {^n-l-k ■ Vx)^fc-lj = 0. (A6) 

Any quantity with a negative subscript is assumed in this paper to be zero. 
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Appendix B. Bounds for a solution to the problem (21) 
for a lineeirly stable CHM state 

Suppose a CHM state V, H, © is linearly stable to small-scale perturbations, which 
is understood here as an exponential decay of any solution to the problem 

M(W) = 0, (Bl) 

Vx-W^^: (w'^)„ = (w'^)^ = (B2) 

in any suitable norm || ■ ||, e.g., in the norm of the space of times differentiable 
functions . Wc show in this appendix, that for such a CHM state any solution 
w = (a;', h', 9) to the problem (21) is globally bounded in this norm in time for any 
permissible initial conditions and right-hand sides f = (f^, f"^, /'^), and, moreover, it 
decays exponentially in fast time, if f does so. (Slow variables play the role of passive 
parameters here, and we do not explicitly refer to them.) 

By the assumption on linear stability of the CHM state, any solution to the 
problem (B1)-(B2) satisfies the inequality 

II w(x, h)\\< Ce-'^(*i-*2) II w(x, h) II (B3) 

for any ti > ^2 > and any initial conditions satisfying (B2), where C and a > are 
constants, independent of ii, t2 and w. 

By linearity, split a solution of (21) into the sum w = W/ -|- w//, where W/ is 
a solution to (B1)-(B2) with the initial condition w/|t=o — w(x, 0), and W// is a 
solution to (21) with w//|t=o = 0. By virtue of (B3), 

||w7(x,t)|| < Ce-^*||w(x,0)||. (B4) 

By the Duhamel principle (see, e.g., Polyanin, 2001) 

= [ w(x, t, t) dr, (B5) 
Jo 

where w(x, t,r) is a solution to (B1)-(B2) for t > r with the initial condition 
w{x,t,T)\t=T = f(x,T). Hence, by (B3), 



||w(x,t,r)||<Ce-'^(*--)||f(x,r)||, 

which implies 

||w77(x,i)||<-max||f(x,T)||. (B6) 

Combination of (B4) and (B6) shows, that w = Wj -|- W// is globally bounded in time 
in the norm || • ||. 

Now, suppose f satisfies the inequality 



|f(x,t)|| <Cfe 



-at 
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where Cf is a time-independent constant (without any loss of generahty we assume 
the same exponent cr > in this inequahty and in (B3) ). Then (B3) imphes 

||w(x,i,T)|| < Ce-^(*-^)Cfe-^^ = CCfC"'^*, 

and thus from (B5) 

||w//(x,t)|| < CCfte-^' < C'e-"'^ 
for any cr' < cr, i.e. ||w77(x, and hence ||w|| decay exponentially. 
Appendix C. A bound for a solenoidal field 

We demonstrate here that if Vx-v = 0, then ( {v(x, '^)}hW=e(xi,x2)) is asymptotically 
smaller than any power of e. 

Let [/] denote in this appendix the average of / over the plane of the fast hori- 
zontal spatial variables Xi,X2- By the definition of {-j/i the horizontal component of 
^-L/2\.{^}h\ vanishes for all X. Averaging the solenoidality condition V^i^jh ~ 
over the fast horizontal variables and integrating it in the vertical direction, find that 
also [({v},,);,] = for all X. Thus it suffices to show that {{{v}h - [{^}h])\:x.=6{xi,x2)) 
is asymptotically smaller than any power of e. 

The equation 

has a mean-free solution, globally bounded together with its derivatives, since the 
average [•] of its left-hand side vanishes. Denote = dcp/dxi. Then by the chain 
rule 



(({v}.-[{v}.])|x=.(.„..)) = ((^ + 



dxi dx2 \dXi 8X2 




'X.=6{xi,X2) ' \ ^ ■'■ X=£(a;i,a;2) 



(Here d/dxk denotes "the complete partial derivative in Xk" .) Since the partial deriva- 
tives have a zero mean over the plane of the fast horizontal variables, the same op- 
eration can be applied to them an arbitrary number of times. This concludes the 
demonstration. 

A similar statement follows from the same arguments: if /(x, t, X, T) is smooth 
and globally bounded with its derivatives, and (/) = 0, then (/|x=£(a;i,i;2)) — 0{e"') 
for any n > 0. 

Appendix D. Exponential decay of the fluctuating pcirt 
of an integral in time. 

Let us show that -^/q ^{t')dt'^ exponentially decays in time in any norm || • ||, in which 
^{t) does. Suppose 

\\m\\<ce-^\ 

where cr > and C are some constants. Then 
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II fU{t')dt'} = f i{t')dt' - \im I r r i{t')dt'dt" 

^ Jo " Jo t->oo t Jo Jo 

= lim i /* /* ^{t')dt'dt" - - lim i /* /* ${t')dt'dt" 

i-^oo t Jo Jt" t-^oc t Jt Jt 



mdt' }} 



ft ft" , fOO 

/ / Ce-'^' dt'dt" < / Ce-^' dt' = {Cla)e-''\ 
Jt Jt Jt 



1 /•* /•* 
< lim ^ / / Ce""* 

t-+oo t Jt Jt 



Appendix E. A simplified form of the combined eddy diffusivity operator 

A simplified expression for the eddy diffusivity operator can be derived for solenoidal 
((voK and ((hoK: 



2 2 2 



VxxEEE 



2 //_ 92 _ g2 _ Q2 



dXl ^''dXt, '^'dX.dX,^ 



dXl 



^ Q2 ^ Q2 \ \ 



dXl_, 



' dXidX2 



where, for k = 1,2 and n = 3 — k, 

DZ = \ {{DZ,.)i - {DZ,2)i - {DZ,i)i + (^2:2,2)2 - {0^2)2 - {DZ,h) , 

-^2,'ifc — {■Dn,k,n)kj -^s'ifc — (-^l,'ifc,2)fe + (-^2,'ifc,l)fe ~ (-^n'n,n)fc ~ {Dn,k,n)n'i 



changing here the first superscript "f" to "/i", obtain the expressions for -0^'^. 
Similarly, 

Vx X E E ^ (d;:;:^^o)). + r>tMho}k) 



m=l k=l 



Vxx 1:1:5^(5 

m=l k=l '-'^rn 



v,h 
m.k 



k + D^m\ihollk 



) efc. 



where, again denoting n — ?> — k, 

Dn!k = {Dn!k)k: D^^k = {Dk,t)k - {Dl'^)k " (-D^',fe)n, 

and expressions for D^'j^ are obtained changing here the first superscript 'V to "/i" . 

Appendix F. Equation for the amplitude of the mean-free neutral mode 

In this appendix one of the main results of the paper is presented, the equation for 
the amplitude of the mean- free neutral mode, cq (see (76)): 
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k=l \ ^^'^ m=l '~>^m <~>^m ' 

k=l \ 

m=l / 



2 2 / Q o 
j=lk=l 

TTi — 1 \ 3 3 I 

+ c^ci + E f (^i;^oK + ciih,),k) + E f CO (czMMkv^)im + ^S. 

fe=l \ m=l \ 

+ C5tH)).((/^oK((/^o)), + CS,((/^oK((/^oK((/^o)),) ) ) 

2 2 2 o3 / 

+ E ( Til- ( 'V^„,„,iC„ + (Af«„,„,jj((i>o»t + N,^„^jih„))t) I 

m=l V^^m \ / 
k=l \ 

2 /)3v7-2 / 2 
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It is convenient to represent its coefficients in tfie terms of tfie 10-dimensional fields 

w, = (o„Vi,H,,e,), 

k — V^k ' ^yfc + 6fc) ^fc j'-'fcJ) ^fc — V^k ' ^fc ' ^A: + ^fe) '-'fc J' 

pT-i),- / ^V,U} I f^V,V f~\V,h y^V,9 \ 

^m,k ~ \y'm,k ^ni, A;, 3^=35 '^m.fc' *^m,fc' ^ m,k ) i 

'Pih,- ff-^v,u) 7^h,v f-^h,h ^h,d \ 'PiC,- f f-^c,u) '^c,v f-^c,h /^c,9\ 

^m,k — \y'm,ki^m,ki^m,ki^m,kj 1 ^k — [^k ^^k ^^k i^k J i 

Y^'" — fv'"''^ Y^'" Y''''* Y^'^^ y'*'' — /'y'^''^ y'*''" y'*''* y^'^\ 

m,k \ m,ki m,ki m,ki m,k J ' m,k \ m,ki m,ki m,ki m,k J ' 

Y*^'' = ( Y*^'"^ Y*^'^ Y'^''^ y^'^ ^ 

m,fej \^ m,fcj' m,fej' m,k,j ) ' 

bilinear forms 

^(p-, q) = (((p'' X q'^ + X p'^ - p'^ X (Vx x q'^) - q" x (Vx x p'^)) • (Vx x S^) 
+ (p^ x q'^ + q^ x p'^) • (Vx X S^J) - ((p'' • V^)q' + (q^ • V^)p') Sl)), 
C,(p-, q) = (((p'' X q'^ + q'' X p'^ - p'^ X (Vx X q'^) - q'^ X (Vx X p^)) X 
+ ((Vx X S:) X p'*) X q'* + (q^ X p'* + p^ x q'^) x - v'q'Sl)), 
of vector fields p' — (p'^, p'', p'*,p^) and q' — (q'^, q^, q'^, q^), and of linear forms 
S{i) = (((f X fio) • (Vx x S:-) + (f X Ho) • (Vx X S^:) - (f • Vx)eo-5f)), 

J-,(p-) = 2 (( Sr+r;^- } + C, (Wo, p) + /3o((A3 x S:')),. 

In this notation, 

= ^,(S-) + £{X;') + i3(Wi, Gf ) + C,(Wi, S-), 

^i:;=-^(sr) +^(^:;) +i3(Wi,G;;,) +c,(Wi,sr), 

= ^,(Q-'-) + £:(A^^) + (C,(S-, S-) + B(S-, Gf )) /2, 

= ^.(QD + + c,(s^.'-,s-) + B(sr,q'-), 

^^'^ = -^.(Qf '■) + si^tn + S-) + ^(s^, Gf ), 

4:/ - .^.(Qf '■) + ^i^tn + C,(S-, S^) + ^(S-, Gj;,), 

^m,k,j — '^J\^m,k ^ ^k,m) ^ '^\-^m,k,j ^ ^k,m,j) ^ '-'jl'^m ^ "^k ) ^ "I'^m ; ^j,k)^ 
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(where is the Kronecker symbol) , 

q'' = B(Wi, S-) + -B(W2, S-) - (52iV^S' X S:))3, 

= B(Wi, S^) + B(W2, ) - MV.Sf X S-))3, 
= i3(Wi,S^) + i3(W2,Si''-) -/32((V.5l''^ x S:))3, 
g- = i3(Wi,Q^^'-) + ^(S-,S-), 
C = i3(Wi, Qr-) + i3(S-, Sr ) + S(S-, sr), 
gf = S(Wi, Qf '■) + B(S-, S^) + S(S-, S^), 

C,.= ^(Wi, Qr;,) +B(s;;;-, s^), ^(Wi, qJ',) + B(sr, st) + B{s^^, s^), 

C = ^(Wi,QS)+^(S:^-,St), 

C"= = -B(Q'^'^'-,S-), 

Q = i3(Qr'-, S-) + S(Q^^'-, Sr), C^, = B{Qt'-, S-) + S(Q^^'-, S^'), 

^m,fe = '^(Qm.'fe' S ) + ^(Qfc ' > ), C^^. = ^(Q^^Ji., S ) + ^(Qfe ' , ), 

Cj, = B(Qj'„ S-) + S(Qr-, S::;-) + B(Q^^'-, S^,''), 

^m,k,i ~ '^V^rn.k^'^j J) ^m,k,j ~ '^V^m,k^'^j )i 

= -^((^r.j)neO + <5^5^(£:(d:';;,) + ^.(n^)), 

^i,n,m,k,j ~ ~^ ii^m,k,j)n^i) + '^^'^01(^(^1, raj) + '^jC^i,n))j 

Nin,m, = -^((-4:;;:;)ne,) + 5^<5)B(Wi, Yr;„), 

^nj = -^((^l)ne.), A7:;„.j = -^((^™:I)rae.), N^Xn., = ^^((^S^)ne.), 
^i,n,m,k,j — ~^ii^m',k,j)n^i), ^ln,m,k,j — ~^ii^m^k,j)n^i), 

^i,n,m,k,j ~ ~^ {{■^m'k,j)n^i) , 
^i,n,m,k,j ~ ~^ {i^m,k,j)n^i) j ^i,n,m,k,j ~ ~ ^{{^m,k,j)n^i) 1 
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^i,n,m,k,j,p ~ ^\\"'m,k,j,p)n^i) j 

i;ni i,m^^ Ji i,m,j i,mi^j li i,m,j i,rn^^j )) 

Mi = s(Yj;;„s-), Mi,^. = s(Yj;;„sj'-), m;^,,=s(y;;„,s,^'-), 

All the coefficients M and N vanish, if the CHM state Vo,Ho,©o possesses a 
symmetry without a time shift {T — 0). 

Appendix G. Evaluation of coefficients of the eddy correction terms 

It suffices to solve 28 auxiliary problems (4, 8, 10 and 6 of types I-IV, respectively) to 
evaluate the coefficients D '', d''' and A "'' in the mean-field equations (59) and (62) 
for the weakly nonlinear stability problem (or these coefficients in (71), (73), (86) and 
(88) ). However, it is possible to exploit the fact that solutions to the auxiliary prob- 
lems (G^_^, G^^^j, G^^/i), {Qm,k^ Qm,fc; Qm,k) 0^m,k^ ^m,h^ ^m,h) enter the aver- 
ages (63) and (64) only in scalar products with the vector field (H x e3, —V x e3, 0), 
and in (60) and (61) with W„j = (-1^-'e„ - Vej^H^e,, + H''ej,0) (4 vector fields 
for n,j = 1,2). We refer to the scalar product ((7 ■ W)) of 7-dimensional vector fields 
7 = (y , and W = (W^ W^, W^) in the layer. 

Since in (59) gradients, on which the curl acts, can be neglected, it suffices to 
evaluate in (60) and (61) scalar products with Wi^2, W2,i and Wi^i — W2,2- The 
number of auxiliary problems to be solved for evaluation of coefficients in (59) and 
(62) is reduced to 8, if auxiliary problems for the adjoint operator are considered, 
following Zheligovsky (2005, 2006a,b). They have the same numerical complexity as 
the auxiliary problems considered above. Then evaluation of coefficients D'^'', d'^'', 
A'\ A and A'^ -' (see (72), (74), (87), (89) ) in new terms in equations (71), (73), 
(86) and (88) also does not require solving auxiliary problems of types V and II'-V: 
it is only necessary to compute the mean-free eigenvector S' e ker A^. 

Let 

Y = {Y)v + r + ^Y, i'' = h'')h + t + ^i^ (Gi) 

(in this appendix all differential operators are in fast variables) be decompositions of 
three-dimensional vector components of W and 7 into spatial mean, solenoidal and 
potential parts. Suppose 7 and Z satisfy the equations 

C%rn^fnn-F\ £^(r,7') = F', C'{j\f)^F', (G2) 

V-r = V-7" = 0, (7''), = (7'^), = 0; 

{C^YiZ", Tl\ Z% = W^ (£^)'XZ", Z'') = W^ (£^)^(Z^ Z^) = O, (G3) 

V ■ Z" = V ■ Z'^ = 0, {Z'")h = {Z'')h = (G4) 

and the boundary conditions for vector fields in the domain of the adjoint operator 
(in our case coinciding with the boundary conditions for the domain of C, (2) and 
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(4) ). Here is the adjoint operator to the restriction of C (defined in section 4) 
onto the subspace of its domain, where horizontal magnetic components have zero 
spatial means: 

£t(v, h, e) = ((r)^(v, h, e), {(£*)'^(v, h, e)}u, (r)^(v, h, e)). 

Then 

((y . + 7^ ■ W^)) (G5) 

= i{^')h ■ {Y)h + {'w'')h ■ {i'')h - w^v ■ Y - w^v ■ 7" 

■ ((>ctr(z^z^z^),(£t)'^(z^z''),(£t)^(z^z^)))) 

= (((W"),, ■ {Y)h + {W^)h ■ {l^)h - ■ Y - w^v . 7^^ + . + • Z'* + F^Z^)) . 

Thus, for evaluation of the eddy coefficients it is enough to solve the auxiliary 
problems for the adjoint operator (G3)-(G4) and to modify the statements of aux- 
iliary problems in order to match the form (G2) by making substitutions (Gl) and 
considering vector potentials of the vorticity equations. Only solutions to auxiliary 
problems of types 11 and 11' are not solenoidal and have non-vanishing spatial means 
of horizontal magnetic components; their divergences are determined by (43e), (44e) 
and (79e) and the spatial means by (56); for the problem (79) the spatial mean can 
be found applying (17). Equations in the form of vector potentials of the vorticity 
equations for problems of type 111 are presented by (53)-(55). They can be also easily 
obtained for problems of all other types, except for types 11 and 11'; for these equa- 
tions, "uncurling" can be performed numerically in the Fourier series representation. 

If the perturbed CHM state V, H, O is steady, then a steady solution to (G3)- 
(G4) is sought. Numerical solution of the evolutionary problem (G3)-(G4) for time- 
dependent data can be problematic, since the operator C'^ ceases to be parabolic. 
However, it becomes parabolic, if time is reversed. Thus, it is possible to overcome 
this difficulty, setting initial conditions for Z '' at t = £ > 0, and solving (G3)-(G4) 
for decreasing t io t — Q. This transformation suffices, if the CHM state V, H, 9 
is periodic or quasiperiodic, and a periodic or quasiperiodic, respectively, solution 
is sought. (Alternatively, such a solution can be found numerically in the form of 
Fourier series in time.) If time dependence is more complex, the spatio-temporal 
averaging of the scalar products involving Z'-' in (G5) can be implemented as 

1 ft i-L/2 ft/2 fi/2 

((F • Z)) = lim lim / / / / F(x,t) ■ Z(t;x,t) dxidx2dx3dt 

i-^ooi-^ootLP Jo J-L/2J-e/2J-e/2 

(if this limit exists), where Z(t;x, t) denotes the solution, obtained with the time 
reversal. 
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